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ABSTRACT 


Two  major  tasks  performed  during  the  report  period  of  investigation  were:  (i^ 
laser  ignition  of  nitramiue-based  composite  propellants  using  a  high-powered  CO2' 
laser  and  (ii)  crack  propagation  and  branching  in  burning  solid  propellants. 

The  laser  ignition  of  a  series  of  RDX-based  composite  propellants  was  studied 
theoretically  and  experimentally.  A  comprehensive  model  for  tho  radiative  igni¬ 
tion  of  the  nitramine  composite  propellants  was  formulated.  The  theoretical  model 
was  solved  numerically.  A  radiative  ignition  test  setup  employing  a  high-powered 
CO?  laser  was  designed  and  constructed  to  study  ignition  characteristics  of  a  series 
of  RDX-based  nitramine  composite  propellants.  Interactions  between  the  infrared 
laser  beam  and  nitramine  composite  propellants  revealed  a  number  of  interesting 
processes.  The  laser  ignition  of  composite  propellants  involves  many  complex  ther- 
mophysiochemical  processes  including  gasification,  initiation  of  a  luminous  flame, 
propagation  of  the  flame,  and  chemical  reactions  near  the  sample  surface  as  well  as 
heat  conduction  in  solid  phases.  In  general,  the  delay  time  based  upon  the  onset 
of  light  emission  decreases  monotonically  with  the  increase  of  incident  laser  energy 
flux.  /The  slope  was  found  to  be  affected  by  the  oxygen  concentration  in  the  ambient 
gas.  This  effect  is  prominent  when  the  sample  contains  volatile  energetic  plasticizers. 
The  eifect  of  oxygen  concentration  is  more  pronounced  in  the  gas-phase  ignition  of 
nitranMne  propellants  than  in  conventional  AP-based  composite  propellants.  The 
flam/ structures  of  nitramine  propellants  depend  strongly  upon  the  ambient  gas 
conditions.  As  ambient  oxygen  concentration  is  increased,  the  flame  moves  closer 
tcx^he  sample  surface;  this  implies  that  the  gas-phase  reaction  with  ambient  oxygen 
lays  an  important  role  in  the  initial  phase  of  laser  ignition  of  nitramine  propellants. 

The  importance  of  crack  propagation  in  solid  rocket  motors  is  widely  recog¬ 
nized.  However,  the  processes  of  crack  propagation  and  branching  in  burning  solid 
propellants  are  not  as  yet  well  understood.  These  processes  could  be  instrumental  in 
creating  large  specific  areas  for  burning,  and  in  extreme  cases  causing  rocket  motor 
failure.  Numerous  samples  have  been  recovered  from  interrupted-burning  experi¬ 
ments.  Depending  on  the  pressurization  rate  and  sample  geometry,  different  modes 
of  crack  propagation  and/or  branching  were  observed.  This  problem  was  analyzed 
through  the  use  of  basic  physical  principles,  current  theories  on  crack  propagation 
in  viscoelastic  media,  and  experimental  results.  A  set  of  governing  dimensionless 
parameters  which  control  and  characterize  the  degree  of  damage  was  obtained.  Key 
factors  influencing  the  crack  combustion,  propagation  and  branching  process  con¬ 
sidered  in  the  dimensional  analysis  include  chamber  pressurization  rate,  initial  crack 
length,  initial  temperature,  sample  geometry,  and  mechanical  and  thermal  proper¬ 
ties  of  the  solid  propellant.  Based  on  experimental  results,  a  multiple  regression 
analysis  was  performed  to  obtain  a  set  of  correlations  between  the  dimensionless 
parameters.  The  procedures  developed  in  this  work  demonstrate  a  methodology  for 
obtaining  similar  correlations  which  could  be  used  as  guidelines  for  safe  operating 
conditions  of  rocket  motors. 


SECURITY  CLASSIFICATION  OF  THIS  PAGEflFhan  Data  Entered) 


LASER  IGNITION 

OF  NITRAMINE  COMPOSITE  PROPELLANTS 

AND 

CRACK  PROPAGATION  AND  BRANCHING 
IN  BURNING  SOLID  PROPELLANTS 


FINAL  REPORT 


Sponsored  by 


Office  of  Naval  Research 
Power  Program 
Arlington,  Virginia 


Contract  No.  N00014-79-C-0762 


Prepared  by 


Accesion  For 

NTIS  CRA&I 
D1IC  TAB 

Unannounced 
Justification 


Kenneth  K.  Kuo 
Jong  U.  Kim 
Jeffrey  Nitnis 
Jesse  Smedly 
Jir-Ming  Char 
Barry  Fetherolf 
Timothy  Snyder 
Claudio  Zanotti 


By . 

Dl--t:ib'.!iic->  / 

Aoai;  ‘ 

Dist 

Prl 


Department  of  Mechanical  Engineering 
The  Pennsylvania  State  University 
University  Pork,  PA  16902 


Reproduction  in  whole  or  in  part  is  permitted 
for  any  purpose  of  the  United  States  Government 


Approved  for  public  release;  distribution  unlimited 


TABLE  OF  CONTENTS 


Page 

ABSTRACT . iii 

TABLE  OF  CONTENT . v 

LIST  OF  FIGURES . viii 

LIST  OF  TABLES . xi 

NOMENCLATURE . xii 

ACKNOWLEDGEMENT . xv 

I.  INTRODUCTION . 1 

II.  LASER  IGNITION  OF  NITRAMINE  COMPOSITE  PROPELLANTS  ....  3 

2.1  INTRODUCTION  . 3 

2.1.1  General  Statement  of  the  Problem . 3 

2.1.2  Previous  Related  Studies . 5 

2. 1.2.1  Radiative  Ignition  -  Experimental . 10 

2. 1.2.2  Radiative  Ignition  -  Theoretical . 12 

2. 1.2.3  Combustion  of  Nitramine  Propellants  . 16 

2.1.3  Scope  of  Present  Work . 16 

2. 1.3.1  Theoretical  Approach . 16 

2. 1.3.2  Experimental  Approach . 16 

2.2  THEORETICAL  APPROACH  . 18 

2.2.1  Description  of  Physical  Model . 18 

2.2.2  Assumptions . 21 

2.2.3  Mathematical  Model  . 22 

2.2.3. 1  Governing  Equations  and  Boundary  Conditions . 22 

2.2. 3.2  Physicochemical  Consideration  for  Pyrolysis 

and  Ignition  of  Nitramine  Composite  Propellants . 26 

2. 2. 3. 3  Modeling  of  Source  Terms  . 31 

2.2. 3. 4  Transformation  of  Governing  Equations . 32 

2.2.4  Numerical  Method . 33 

2.3  EXPERIMENTAL  APPROACH . 34 

2.3.1  Introduction . 34 

2.3.2  Experimental  Apparatus  and  Operation . 34 


2.3.3  Data  Acquisition  and  Analysis  Systems . 40 

2.3.4  Propellant  Samples . 41 

2.3.5  Beam  Intensity  Cistribution . 41 

2.3.6  Temperature  Measurement . 45 

2.3.7  Gas  Species  Measurement . . . 49 

2.3.8  Surface  Structure  Observation . 49 

2.3.9  Schlieren  Technique . 50 

2.4  RESULTS  AND  DISCUSSION  . 52 

2.4.1  Experimental  Results . 52 

2.4.1.1  Ignition  Delay  Time . 52 

2.4. 1.2  Flame  Observation . 58 

2.4. 1.3  Thermal  Wave  Structure . 75 

2.4. 1.4  Gas  Analysis . 81 

2. 4. 1.5  Sample  Surface  Structure . 84 

2.5  SUMMARY  AND  CONCLUSIONS . 87 

CRACK  PROPAGATION  AND  BRANCHING  IN  BURNING  SOLID 
PROPELLANTS  . 89 

3.1  INTRODUCTION  . 89 

3.1.1  Motivation  and  Objectives  . 89 

3.1.2  Brief  Survey  of  Related  Work . 90 

3.2  EXPERIMENTAL  SETUP . 93 

3.2.1  Selection  of  Sample  Propellant  Type . 93 

3.2.2  Sample  Configuration . 94 

3.2.3  Test  Setup . 94 

3.2.3. 1  Test  Chamber . 94 

3.2.3.2  Data  Acquisition  System .  99 

3.2.4  Test  Procedure . 102 

3.3  THEORETICAL  ANALYSIS . 103 

3.3.1  Crack  Sample  Geometry . 103 

3.3.2  Determination  of  Mechanical  Properties . 103 

3.3.2. 1  Relaxation  Modulus . 105 

3. 3. 2. 2  Crack  Propagation  Velocity . 105 

3. 3. 2. 3  Fracture  Energy  and  Related  Parameters  . 110 

3.3.3  Fnergy  Balance  Applied  to  the  Control  Volume . Ill 

3.3.4  Dimensional  Analysis  . . . . . 114 


vii 


I 

8 

I 

1 

I 

1 

8 


3.3.4. 1  Arrangement  into  Meaningful  Dimensionless 


Groups . 114 

3.3.4.2  Functional  Relationships  Between 

Dimensionless  Parameters . 115 

3.4  RESULTS . 117 

3.4.1  Determination  of  Dimensionless  Parameters . 117 

3.4. 1.1  Pressure  Loading  in  the  Damaged  Zone . 117 

3.4. 1.2  Interpretation  of  Film  Data . 121 

3.4. 1.3  Determining  Other  Key  Parameters . 121 

3.4.2  Relationships  Between  Dimensionless  Parameters . 127 

3.4.3  Discussion  of  Results . 129 

3.5  SUMMARY  AND  CONCLUSIONS . 135 


REFERENCES 


137 


1 

i 

i 


APPENDICES: 


A.  MEASUREMENT  OF  C02  LASER  ENERGY  FLUXES  . . 144 

B.  DIFFERENT  TECHNIQUES  USED  FOR  ACHIEVING 

BEAM  UNIFORMITY . 149 

C.  ABSORPTION  OF  C02  LASER  BEAM 

ENERGY  BY  C02  GAS . 154 

D.  OPERATION  PROCEDURE  AND  SAFETY  FEATURES  OF 

THE  C02  LASER  SYSTEM . 158 

E.  SPECIFICATION  OF  THE  EQUIPMENT  USED  IN  THE 

IGNITION  STUDY  . 163 

F.  TEST  PROCEDURES  AND  CHECKLIST  FOR  CRACK 

PROPAGATION . 167 

G.  ERROR  ANALYSIS . 172 


n 
& 


viii 


Figure 

1. 

2. 

3. 

4. 


5. 

6. 

7. 

8. 

9. 

10. 
11. 
12. 

13. 

14. 

15. 

16. 

17. 

18. 

19. 

20. 
21. 
22. 
23. 


LIST  OF  FIGURES 

Page 

Physical  and  Chemical  Processes  in  Composite  Solid  Propellant 


Ignition  under  Radiative  Heating  . 19 

Heterogeneous  Structure  of  Nitramine  Composite  Propellants  . 20 

Statistically  Averaged  Element  and  Finite  Difference 

Grid  Pattern  . 20 

Physical  and  Chemical  Processes  in  the  Ignition 

and  Combustion  of  Nitramine  Composite  Propellants  . 27 

A  Schematic  Diagram  of  COa  Laser  Ignition  Test  Setup  . 35 

Low  Pressure  Chamber  . 37 

A  Schematic  Diagram  of  Gas  Control  System  . . . 38 

The  CO2  Laser  . 39 

Configuration  of  Fine  Wire  Thermocouple  Mounting . 46 

The  Junction  of  Thermocouple  (Pt-Ptl3%Rh,  25 pm)  . 47 

Schematic  Diagram  of  Temperature  Measurement  . 48 

Schlieren  Optical  System  . 51 

Delay  Time  for  Onset  of  Light  Emission  (BLX4)  . 54 

Delay  Time  for  Onset  of  Light  Emission  (BLX8)  . 55 

Delay  Time  for  Onset  of  Light  Emission  (BLX9)  . 56 

Delay  Time  for  Onset  of  Light  Emission  (AP/HTPB)  . 57 

Ignition  Characteristics  of  BLX8(59603,lartm)  . 58 

Flame  Development  of  BLX8  ( q "  —  7OOW/em2,5%02)  . 61 

Flame  Development  of  BLX4  ( q "  =  700 W/cm2, 1  %Oi)  . 62 

Flame  Development  of  BLX4  (y"  =  700W/cm2,21%03)  . 63 


Flame  Structure  Recorded  by  Conventional  Video  Camera . 65 

Schlieren  Pictures  and  Direct  Pictures  . . .  66-74 

Surface  Temperature  vs.  Time  (q"  =  1000 W j cm2,  He)  . 76 


ix 


25.  Surface  and  Gas  Temperature  vs.  Time  (qn  =  lOOOW/em2)  . 78 

28.  Surface  and  Gas  Temperature  vs.  Time  (gM  =  8001V/cma)  . 79 

27.  Gas  Temperature  vs.  Time  ( q "  =  800W/em2)  . 80 

28.  A  Typical  Output  of  Gas  Chromatograph 

for  Porapak  Q  and  R  Column . 82 

29.  A  Typical  Output  of  Gas  Chromatograph  for  MS  5A  Column  . 83 

30.  Scanning  Electron  Micrographs  of  BLX4  . 86 

31.  Crack  Propagation  Velocity  versus  Pressurization  Rate . 91 

32.  Propellant  Sample  Configuration  Before  Loading  . 95 

33.  Test  Chamber  . 96 

34.  Driving  Motor  . 98 

35.  Data  Acquisition  System  . 100 

36.  Geometry  of  the  Damaged  Propellant  Sample . 104 

37.  Relaxation  Modulus  Time  Shift  Factor  versus  Temperature  . 106 

38.  Master  Relaxation  Curve  . 107 

39.  Test  Configuration  for  Crack  Velocity  Characterization . 108 

40.  Crack  Propagation  Velocity  versus  Stress  Intensity 

Factor  for  ARC  4525  ...  .  109 

41.  Control  Volume  Used  in  Energy  Analysis  . 112 

42.  Pressure-Time  Trace  for  DNICP  27  . 118 

43.  Recovered  Sample  from  DNICP  27  . 119 

44.  Pressure-Time  Trace  for  DNICP  40  . 120 

45.  Film  Record  from  DNICP  40  . 122 

46.  Instantaneous  Bright  Film  Contour  From  DNICP  40  . . 124 

47.  Crack  Tip  Location  versus  Time  for  DNICP  40  . 125 

48.  Instantaneous  Bright  Film  Contour  From  DNICP  31  . 124 

49.  Data  Correlation  for  the  Extended  Crack  Length  Tlx  . 130 

50.  Data  Correlation  for  the  Rate  of  Increase  in  Specific 

Surface  Area  H©  . 131 


Data.  Correlation  for  the  Number  of  Macrocracks  Tic  •  • 
Data  Correlation  for  the  Crack  Propagation  Velocity  flo 


xi 


LIST  OF  TABLES 


Table  Page 

I.  Composition  of  Test  Propellant?  . 42 

II.  Thermophysical  Properties  of  the  Ingredients  . 43 

III.  Molecular  Structure  of  Major  Ingredients . 44 

IV.  Film  Interpretation  of  Crack  Propagation  Processes  . 123 

V.  Steady-State  Burning  Rate  Constants . 127 

VI.  Experimental  Value  of  Dimensionless  Parameters  . 128 

VII.  Estimation  of  measurement  Errors . 173 


MM 


""-$5 


xii 


NOMENCLATURE 


Symbol 

Afc 

An 

Aa 

a 

<»o 

or 

6 

e,ep 

P 

E 

9 

hp 

At 

/ 

k 

ki 

K, 

Kin 

L 

Ld 

MW 

m 

rhft 

ms 

n 

P 


Description 

Burning  surface  area  in  the  affected  region 

Total  burning  surface  area 

Specific  burning  surface  area 

Total  crack  length 

Initial  crack  length 

Time  shift  factor 

Propellant  sample  thickness 

specific  heat 

binary  mass  diffusion  coefficient 

activation  energy,  relaxation  modulus  or  total  energy 

acceleration  of  gravity 

Enthalpy  of  the  propellant 

Stagnation  enthalpy  of  product  gases 

radiation  flux 

thermal  conductivity 

rate  constants  for  reaction  i 

Stress  Intensity  Factor 

Critical  stress  intensity  factor 

thickness  of  cylindrical  pellet 

Length  of  extended  macrocrack 

mean  molecular  weight  of  gas  mixture 

Mass  flow  rate 

mass  burning  rate  of  RDX  filler 
mass  burning  rate  of  binder 
Number  of  macroscracks 
pressure 


v 


V  V 


Q,pyo 

Q 

r 

Tb 

R 

Ru 

Rl,R2 

J23 

T 

t 

VM 

V 

Va 

Vc 

Vp 

Yf 

z 

Z 


heat  flux  (energy  per  unit  time  per  unit  area) 
heat  generation  rate  (energy  per  unit  time  per  unit  volume) 
heat  of  pyrolysis  per  unit  mass 
Heat  loss 

radiat  distance  from  center  of  pellet 
burning  rate 
gas  constant 
universal  gas  constant 

radii  of  outer  surfaces  of  filler  and  binder,  respectively 
Sample  multiple  coefficient  of  determination 
temperature 
time 

gas  phase  velocity  in  z  direction 
Voltage 

Affected  volume 
Chamber  free  volume 
Crack  propagation  velocity 
mass  fraction  of  species  j 

distance  from  gas-solid  interface;  positive  in  gas 
pre-exponential  factor 


Greek  Letters 

a  absorptivity  or  thermal  diffusivity 

/?  in-depth  radiation  absorption  coefficient 

7  Specific  heat  ratio 

T  Fracture  energy 

c  emissivity  or  strain 

v  Poisson’s  ratio 

£  Shifted  time 


p 

a 


& fm 
T 


V>o 


density 

Stress 

Maximum  stress  in  the  failure  zone 
transmissivity  of  solid  phase 
Crack  cavity  angle 
Outside  taper  angle 
mass  production  rate  per  unit  area 
mass  production  rate  per  unit  volume 


Subscripts 


1,2, 3, 4, 5 

B 


9 

i 

pyro 


various  reactions  of  species 

binder 

gas  phase 

initial  value  or  species 
pyrolysis 


R  RDX  filler 


rad  radiation 

s  solid  phase 

solid-gas  interface 


9-9 


XV 


ACKNOWLEDGEMENT 


This  research  has  been  sponsored  by  the  Office  of  Naval  Research,  Arlington, 
Virginia,  under  Contract  No.  N00014-79-C-0762.  Dr.  Richard  S.  Miller  served  as 
the  technical  monitor  and  program  manager  for  this  contract.  His  support  of  this 
research  investigation  is  greatly  appreciated.  We  would  also  like  to  acknowledge  the 
encouragement  of  Mr.  Tom  Boggs  of  NWC  for  investigating  nitraroine  propellant 
ignition.  Appreciation  is  extended  to  Drs.  E.  C.  Martin  and  Rena  Y.  Yee  of  the 
Naval  Weapons  Center  for  providing  various  types  of  nitramine  propellants  for  this 
ignition  study.  Thanks  are  also  due  to  Professor  R.  A.  Schapery  of  Texas  A.  &:  M. 
University  for  usefuj  exchange  of  ideas  in  crack  propagation  studies.  The  authors 
would  also  like  to  thank  Mr.  T.  Torikai  of  Ashahi  Chemical  Ind.  Co.  of  Japan  for 
his  assistance  on  this  study. 


I.  INTRODUCTION 


This  is  the  final  report  summarizing  the  research  investigation  conducted  dur¬ 
ing  January  1,  1986  to  August  31,  1987,  under  the  project  entitled  “Ignition  of 
Solid  Propellants  and  Propagation  of  Burning  Propellants  Cracks”  (Contract  No. 
N00014-79-C-0762) . 

The  major  objectives  of  this  investigation  are:  1)  to  develop  a  comprehensive 
theory  for  the  radiative  ignition  of  nitramine  composite  propellants  and  to  evaluate 
the  theory  by  comparison  with  experimental  results  obtained  in  this  investigation 
and  2)  to  examine  propagation  and  branching  processes  in  burning  solid  propellant 
crack  samples.  Specific  objectives  of  this  study  are: 

1)  to  determine  the  dependence  of  delay  time  for  onset  of  light  emission  (tLE)  as 
a  function  of  incident  energy  flux  (<?"),  ambient  gas  compositions,  and  ambient 
pressure; 

2)  to  observe  the  site  of  initial  light  emission  and  flame  development  processes; 

3)  to  measure  gas-phase  temperatures  in  the  vicinity  of  the  burning  surface; 

4)  to  study  the  surface  structures  of  test  samples  recovered  before  and  after  onset 
of  light  emission; 

5)  to  compare  the  ignition  characteristics  of  a  series  of  nitramine  solid  propellants 
with  those  of  a  conventional  AP/HTPB  propellant; 

6)  to  develop  a  theoretical  model  and  to  obtain  numerical  solutions  for  ignition 
of  nitramine  composite  propellants  under  CO2  laser  heating; 

7)  to  validate  the  theoretical  model  with  experimental  data; 

8)  to  study  important  parameters  governing  the  laser  ignition  of  nitramine  com¬ 
posite  propellants; 
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9)  to  define  a  group  of  meaningful  dimensionless  parameters  for  characterizing  the 
operating  conditions  and  resulting  damage  of  a  solid  sample; 

10)  to  obtain  test  data  necessary  for  developing  correlations  between  the  dimen¬ 
sionless  parameters;  and, 

11)  to  establish  correlations,  demonstrating  a  methodology  for  generating  useful 
safety  criteria. 


II.  LASER  IGNITION  OF  NITRAMINE  COMPOSITE  PROPELLANTS 


2.1  INTRODUCTION 

2.1.1  General  Statement  of  the  Problem 

Recently,  solid  propellant  researchers  have  focused  a  good  deal  of  attention 
on  the  study  of  nitramine  composite  propellants.  Nitramine  composite  propellants 
offer  several  advantages:  (1)  high  energy,  (2)  low  combustion  temperature,  (3)  low 
infrared  emissions  in  the  plume  region,  and  (4)  low  vulnerability.  However,  most 
studies  [1-11]  have  been  focussed  on  steady-state  combustion.  Little  effort  has  been 
directed  to  the  study  of  ignition  of  nitramine  composite  propellants[12,13]. 

Radiative  ignition  techniques  have  been  used  widely  for  research  and  for  test¬ 
ing  the  ignition  characteristics  of  newly  formulated  propellants[14].  The  radiative 
energy  sources  are  usually  produced  from  an  arc  image  furnace  or  a  COj  laser.  The 
arc  image  furnace  was  a  popular  device  for  laboratory  radiative  ignition  studies 
prior  to  development  of  the  CO3  laser.  In  recent  years,  the  COj  laser  has  been 
more  widely  used  in  propellant  ignition  studies[l5-18,28,29]  due  to  several  major 
advantages  it  offers  over  the  arc  image  furnace[9].  The  advantages  of  employing  a 
CO3  laser  system  for  the  study  of  pyrolysis  and  ignition  are:  (1)  the  energy  flux 
to  the  propellant  sample  can  be  controlled  with  a  high  degree  of  reproducibility; 
(2)  the  environmental  parameters  ,  such  as  chamber  pressure,  initial  temperature, 
and  chemical  composition  of  surrounding  gases  in  a  CO3  laser  ignition  device,  can 
be  varied  independently;  and  (3)  the  improved  CO2  laser  beam  intensity  can  reach 
a  sufficiently  high  level  to  simulate  the  broad  range  of  energy  fluxes  in  propulsive 
devices. 


The  low  vulnerable .  characteristics  of  composite  propellants  used  in  an  am¬ 
munition  system  are  important  factors  in  screening  newly  formulated  propellants. 
Low  susceptibility  to  external  thermal  energy  stimulation  is  a  desirable  feature  of  a 
Low  Vulnerability  Ammunition  (LOVA)  propellant.  LOVA  propellants  should  also 
possess  proper  ignitibility  for  propulsive  purposes.  Recently,  Roller  et  al.[13]  found 
that  LOVA  propellants  can  be  ignited  more  effectively  using  oxidizer-rich  ignition 
materials.  This  further  substantiates  the  importance  of  the  effect  of  ambient  oxi¬ 
dizer  concentration  in  the  ignition  processes  of  LOVA  propellants.  These  ignition 
characteristics  can  be  determined  readily  and  accurately  by  using  a  high-powered 
COa  laser  with  well-controlled  ambient  gas  composition. 

Although  radiative  ignition  of  solid  propellants  drew  extensive  attention  for 
more  than  three  decades,  it  has  its  own  inherent  drawbacks.  Ohlemiller  and  Sum- 
merfield  [14]  critically  reviewed  the  radiative  ignition  of  solid  fuels.  They  found  the 
most  prominent  drawback  to  be  the  lack  of  hot  gas  adjacent  to  the  heated  surface; 
this  results  in  loss  of  the  energy  absorbed  by  the  solid  surface  through  conductive 
heat  transfer  to  the  adjacent  cool  ambient  gas.  Another  major  problem  is  that 
the  radiative  energy  absorbed  by  the  solid  propellants  is  strongly  dependent  on  the 
optical  properties  of  propellants.  Despite  these  disadvantages,  radiative  ignition  is 
attractive. 

The  present  study  emphasizes  achieving  a  better  understanding  of  thermal 
pyrolysis,  flame  initiation,  and  flame  development  in  the  ignition  of  nitramine  com¬ 
posite  propellants. 
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An  extensive  review  of  experimental  and  theoretical  studies  related  to  the  ra¬ 
diative  ignition  of  nitramine  composite  propellants  is  given  in  this  section.  Price  et 
al.[l9]  reviewed  the  state  of  the  art  in  solid  propellant  ignition  theories  published  up 
to  1966.  They  also  discussed  experimental  works  related  to  theoretical  approaches. 
Kulkarni  et  al.(20]  also  reviewed  theoretical  and  experimental  works  on  the  ignition 
of  solid  propellants.  Recently,  Hermance[2l]  presented  an  excellent  review  covering 
work  done  after  Price’s  1966  review. 

Experimental  studies  of  radiative  ignition  will  be  reviewed  first  in  the  present 
study,  followed  by  an  examination  of  previous  theoretical  approaches  to  radiative 
ignition.  Due  to  the  importance  of  information  from  steady-state  combustion  of 
nitramine  propellants  to  the  ignition  study  of  nitramine  composite  propellants,  both 
experimental  and  theoretical  studies  on  the  combustion  of  nitramine  propellants  will 
be  reviewed.  Finally,  since  thermal  decomposition  of  nitramines  is  an  important 
process  in  the  ignition  of  nitramine  composite  propellants,  a  brief  review  of  thermal 
decomposition  of  nitramines  aud  binders  will  be  given  in  this  section  and  will  be 
discussed  further  in  a  later  section. 

2.1.2.1  Radiative  Ignition  -  Experimental 

Shannon[22]  studied  ignition  characteristics  of  composite  solid  propellants  us¬ 
ing  an  arc-imaging  furnace  as  a  radiant  energy  source.  He  used  AP-based  composite 
propellants  with  different  binders  as  test  samples.  A  go/no-go  test  was  used  as  an 
ignition  criterion.  The  pressure  range  was  0.1-4  atm  and  the  range  of  heat  fluxes 
was  10  -  100  cal/cm3  s.  He  defined  the  ignition  process  by  considering  two  sepa¬ 
rate  stages:  (1)  propellant  inert  heating,  and  (2)  ignition  of  the  volatile  following 


propellant  heating.  Thus,  the  ignition  event  can  be  divided  into  time  intervals 
representing  the  two  stages, 

=  tr  +  tc  (2  —  1) 

where 

ti?n  =  total  time  to  ignition 

tr  —  time  of  propellant  heating  by  external  energy  source  alone 

te  time  of  heat  addition  by  various  chemical  reactions  plus  time  of  transition 
to  steady-state  combustion. 

The  rate  of  volatile  generation  rises  rapidly  with  temperature  and  will  exceed 
any  threshold  valu*  within  a  small  temperature  change.  Therefore,  ignition  will 
usually  occur  when  the  solid  reaches  some  fixed  temperature,  and  the  temperature 
will  be  substantially  independent  of  heating  rate  and  other  parameters  over  a  lim¬ 
ited  range.  Shannon  suggested  that  if  a  complete  understanding  of  the  ignition 
mechanism  is  to  be  attained,  the  kinetic  processes  affecting  the  chemical  time,  tc, 

rs 

must  be  defined  in  greater  detail.  The  purpose  of  his  study  was  to  explore  the  inter¬ 
relationship  among  ignition  characteristics  and  propellant  compositional  factors  in 
order  to  determine  chemical  reaction  sequences  that  may  be  important  in  determin¬ 
ing  tc.  He  determined  the  minimum  exposure  times  by  using  the  go/no-go  criterion. 
At  low  pressures,  where  ignition  requires  long  exposure  times,  the  exposure  time 
was  held  constant  and  go/no-go  limits  of  pressures  were  also  determined. 

Shannon[22]  found  that  interfacial  oxidation  between  oxidizer  and  binder  was 
not  important,  and  that  gas-phase  compositional  and  pyrolysis  differences  were 
more  important  than  solid-phase  or  surface  reactions.  He  concluded  that  gas-phase 
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reaction  establishment  is  an  important  process  in  composite  propellant  ignition  by 
radiant  energy. 

Summerfield  and  his  &ssociat^s(l4-18]  conducted  a  series  of  experiments  to 
understand  the  behavior  of  composite  propellants  under  radiative  heating.  They 
suggested  an  equation  for  prediction  of  ignition  delay  time  in  radiative  ignition 
of  composite  propellants.  They  considered  three  dominant  factors  (conduction, 
in-depth  radiation,  and  diffusion/reaction)  in  ignition  delay,  based  on  the  value 
of  radiant  flux.  At  high  energy  flux,  the  delay  time  for  chemical  and  diffusion 
processes  in  the  gas  phase  is  dominant.  At  low  energy  flux,  condensed-phase  heating 
dominates  the  ignition  delay.  In  the  upper  range  of  the  low  energy  flux,  the  in-depth 
absorption  of  radiation  is  more  important  t.ian  heat  conduction.  From  the  above 
reasoning,  they  proposed  a  simple  equation  for  prediction  of  ignition  delay  time. 

tign  =  tco  +  ce(pCAT7(l  -  r)F\  +  (7r/r/4)[pCAr/(l  -  r)F]2  (2  -  2) 

where  a  is  the  sample  absorption  coefficient  for  the  incident  radiation,  r  the  surface 
reflectance,  K  the  thermal  diffusivity,  F  the  incident  flux,  and  AT*  the  temperature 
rise  of  the  surface  necessary  to  yield  ignition.  They  did  not  elaborate  on  any  method 
to  obtain  tcD  *n  Eq.  (2-2).  Thus,  the  equation  is  only  qualitatively  useful  when 
processes  subsequent  to  heating  of  the  condensed  phases  dominate  the  ignition  delay. 

Ohiemiiler  and  Summerfield[15]  examined  the  radiative  ignition  behavior  of 
polymeric  fuels  with  speciU  reference  to  solid  propellants.  Polystyrene  and  an 
epoxy  polymer  in  oxyaen/nitrogen  mixtures  were  studied.  Ignition  delay  time  was 
measured  as  a  function  of  radiant  flux,  pressure,  oxygen  percentage,  and  fuel  ab¬ 
sorptivity.  A  continuous  wave  CO2  laser  was  employed  as  an  energy  source  to 
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eliminate  some  of  complexities  related  to  an  arc  image  furnace.  They  found  that 
as  the  oxygen  partial  pressure  in  the  ambient  gas  increases,  the  ignition  delay  time 
decreases.  At  the  upper  limit  of  oxygen  partial  pressure,  ignition  delay  becomes  in* 
dependent  of  this  variable.  As  the  oxygen  partial  pressure  is  decreased,  the  ignition 
delay  time  increases  until  it  approaches  a  no-ignition  limit.  They  also  studied  the 
effect  of  a  surface  opacifier(carbon  black)  on  ignition  delay.  They  found  that  the 
addition  of  the  opacifier  not  only  reduced  ignition  delay  time,  but  made  the  slopes 
of  ignition  delay  vs.  incident  energy  graphs  steeper. 

Ohlemiller  et  al.(l6]  extended  their  radiative  ignition  studies  to  investigate  the 
dynamic  response  of  a  solid  propellant  to  rapidly  varying  radiative  fluxes.  They 
found  that  when  a  CO3  laser  is  used  as  an  energy  source,  a  flame  was  sustained 
only  during  the  external  laser  heating,  and  was  extinguished  when  the  incident  laser 
beam  is  cut  off.  They  also  found  that  the  tendency  to  extinction  upon  deradiation 
is  reduced  by  increasing  pressure  or  increasing  the  time  interval  for  reducing  the 
energy  flux  to  sero.  They  proposed  a  mechanism  to  explain  this  phenomenon.  A 
short  deradiation  interval  produces  a  corresponding  rapid  decrease  in  surface  gasifi¬ 
cation  rate  and  a  relatively  fast  decrease  in  heat  feedback  from  the  flame.  However, 
the  condensed  phase  acts  as  if  it  has  inertia,  and  cannot  respond  instantaneously 
to  a  change  in  heat  feedback  from  the  flame.  Accordingly,  the  high  condensed- 
phase  temperature  gradient,  previously  established  by  the  external  radiation,  per¬ 
sists  and  causes  a  relatively  large  fraction  of  heat  to  be  conducted  out  of  the  gas 
phase.  This  diminishes  the  flame  intensity.  Thus,  there  is  a  competition  between 
condensed-phase  conduction  processes  that  decrease  the  temperature  gradients  at 
the  condensed-phase  surface,  and  the  weakened  ability  of  the  flame  to  be  sustained 
as  the  condensed  phase  adjusts  to  the  reduced  heating  rate.  They  also  proposed  a 
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theoretical  model  to  simulate  the  extinction  of  the  flame  with  removal  of  Incident 
laser  energy. 

In  the  two-part  study,  DeLuca  et  al.[17-18]  examined  the  ignition  response  of 
several  propellants  under  radiative  heating  and  the  influence  of  the  radiation  source 
(arc  image  vs.  laser)  on  observed  ignition  behavior.  In  the  first  part  of  the  study[17], 
they  emphasised  the  influence  of  propellant  formulation  on  ignition  behavior.  They 
found  that  comparisons  of  thy  chemical  factors  in  the  formulation  can  only  be  made 
properly  when  the  optical  factors  are  minimised.  They  also  found  that  the  HMX 
composite  was  most  difficult  to  ignite  among  the  propellant  samples  tested,  and 
followed  by  AP-based  composite  and  double-base  propellants.  In  the  second  part  of 
the  study[18],  pre-ignition  events  and  effects  of  the  nature  of  the  radiation  source 
were  emphasised.  They  found  that  fast  deradiation  extinction  occurs  only  with  the 
laser,  not  with  the  arc  image  furnace. 

Kashiwagi[23-25]  used  a  COa  laser  in  a  series  of  experiments  to  study  radia¬ 
tive  ignition  of  polymethylmethacrylate  (PMMA)  and  red  oak.  He  found  strong 
attenuation  of  incident  laser  radiation  by  the  plume  consisting  of  decomposition 
products  in  the  gas  phase.  He  postulated  that  PMMA  ignites  by  absorption  of 
incident  radiation  by  decomposition  products  in  the  gas  phase.  He  also  found  that 
the  attenuation  of  radiation  caused  by  decomposition  products  in  the  gas  phase  was 
significant  enough  to  affect  surface  temperature. 

Kashiwagi  and  Ohleroiller[26]  studied  oxygen  effects  on  nonfiaming  transient 
gasification  of  PMMA  and  PE  during  thermal  radiation.  The  presence  of  oxygen 
in  the  gas  phase  was  found  to  increase  the  mass  flux  from  PMMA  and  PE  signifi¬ 
cantly.  However,  the  effect  of  oxygen  on  mass  flux  was  apparently  reduced  as  the 
counterflow  of  decomposition  gases  grew.  Kashiwagi  and  Ohlemiller[27]  also  used  a 


CO)  laser  to  study  the  ignition  mechanism  of  a  liquid  fuel  under  radiation  heating. 

Mutoh  et  al.(30]  conducted  an  experimental  Jtudy  to  explore  the  behavior  of 
the  radiative  ignition  of  PMMA.  A  COj  laser  was  used  as  an  energy  source.  They 
found  a  critical  energy  flux  above  which  ignition  occured  at  the  piume  axis  away 
from  the  PMMA  surface,  and  below  which  it  occurred  at  the  edge  of  the  plume  near 
the  sample  surface.  They  suggested  that  determining  a  priori  the  ignition  criterion 
is  impossible,  and  that  it  is  necessary  to  explore  the  heat  and  mass  fluxes  throughout 
the  solid  and  gas  phases  in  order  to  predict  ignition.  Beckel  and  Matthews[33]  also 
used  a  CO)  laser  to  study  the  radiative  ignition  of  polyxymethylene  (POM).  They 
found  the  effect  of  oxidiier  supply  in  the  ambient  gas  on  ignition  delay  time  of  POM 
to  be  significant. 

Recently,  the  application  of  high  power  lasers  to  gun  systems  as  practical  initia¬ 
tion  devices  has  been  attempted [3 1,32).  Laser  ignition  systems  for  large  caliber  guns 
were  found  to  offer  possibilities  for  improved  performance  and  reliability.  Ostrowski 
et  ai.[31-£?)  employed  a  laser  to  study  ignition  characteristics  of  gun  propellants. 
They  used  a  pulsed  ruby  laser  to  ignite  the  propellants.  They  found  that  proper 
adjustment  of  the  laser  energy  deposition  rates  produced  a  certain  degree  of  self¬ 
confinement  of  the  flame  which  then  allowed  sustained  combustion  to  occur. 

2.1.2.2  Radiative  Ignition  -  Theoretical 

General  ignition  theories  have  been  reviewed  extensively  by  Price  at  al.[l9]  and 
Hermance[21].  In  this  section,  some  recent  radiative  ignition  models  are  reviewed. 

Kashiwagi(34]  proposed  a  theoretical  model  describing  radiative  ignition  of  a  . 
solid  fuel.  He  included  the  effects  of  gas-phase  reaction  and  in-depth  absorption  of 
the  incident  radiation  by  the  solid  fuel.  From  the  results  of  his  theoretical  inves- 
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tic  At  Ion,  he  established  a  range  of  pyrolysis  activation  energies  for  ignition  as  well 
as  the  importance  of  gas-phase  reaction  and  in-depth  radiation.  He  also  found  that 
the  effects  of  ambient  inert  gas  on  the  ignition  delay  time  are  significant.  His  model 
demonstrates  some  important  aspects  of  radiative  ignition;  however,  the  application 
is  limited  to  pure  fuel  ignition. 

Kumar  and  Hermance[35]  presented  a  radiative  Ignition  model  of  a  homoge¬ 
neous  propellant.  The  decomposed  gas  from  the  solid  phase  is  a  homogeneous 
mixture  of  fuel  and  oxidiser.  They  found  that  the  gas-phase  reactions  are  impor¬ 
tant  in  establishing  a  gas-phase  flame,  but  their  influence  in  evaluating  the  ignition 
behavior  of  the  propellant  is  negligible.  They  suggested  that  only  at  very  low 
pressures  and  extremely  high  heating  rates  do  the  gas-phase  processes  become  im¬ 
portant.  Although  this  model  can  handle  the  ignition  of  a  homogeneous  propellant, 
the  applicability  to  a  composite  propellant  is  questionable. 

Kumar[36]  further  extended  the  radiative  ignition  model[35j  to  include  the  two- 
dimensional  nature  of  composite  propellants.  In  this  model,  the  diffusion  of  species 
and  energy  was  considered  two-dimensional.  The  propellant  was  considered  to  be 
opaque.  Thus,  no  in-depth  radiative  heating  was  considered.  Although  the  predic¬ 
tion  of  this  model  agreed  with  some  of  the  experimental  observations  in  other  cases, 
predictions  regarding  the  pressure  dependence  of  ignition  delay  in  inert  atmospheres 
did  not  agree. 

Phuoc  and  Durbetaki(37]  formulated  a  set  of  boundary  layer  equations  which 
included  the  absorption  of  external  radiation  by  the  gas  phase.  They  used  the 
empirical  equations  derived  from  Kashiwagi’s  experiments[24).  They  couid  predict 
the  discontinuity  in  ignition  time  observed  in  Mutoh’s  experiments[30j. 

Kindelan  and  Williams(38]  employed  an  asymptotic  method  for  ignition  of  a 


solid  fuel.  They  included  in  their  analysis  the  endothermic  gasification  of  the  solid 
fuel  and  exothermic  reactions  in  the  gas  phase.  In-depth  absorption  of  the  external 
radiant  energy  flux  was  also  considered.  The  advantage  of  using  an  asymptotic 
method  is  the  availability  of  parametric  studies.  However,  the  applicability  of  the 
asymptotic  technique  is  limited,  especially  when  the  properties  of  the  materials  are 
dependent  on  the  other  variables. 

2. 1.2.3  Combustion  of  Nitramine  Propellants 

McCarty  et  al.[l]  conducted  a  study  to  determine  the  combustion  behavior  of 
nitramine  propellants.  They  studied  the  effects  of  solid  loading  and  particle  size 
on  the  burning  rate  of  HMX  composite  propellants.  They  found  a  slope  break  on 
the  burn-rate  curve  near  136  atm.  They  also  found  that  with  coarse  particles,  the 
HMX  protruded  from  the  quenched  surface,  while  fine  particles  produced  a  smooth 
surface.  This  was  attributed  to  ignition  delay  of  coarse  particles.  With  coarse 
particles,  the  heat  transfer  between  particles  is  reduced  due  to  the  relatively  longer 
distance  between  particles. 

Kubota  et  al.[2]  studied  the  combustion  characteristics  of  nitramine  and  AP 
composite  propellants.  They  found  that  the  replacement  of  AP  with  nitramines 
decreased  the  burning  rate.  They  also  observed  that  the  effect  of  the  particle 
size  of  nitramines  on  burning  rate  was  not  significant.  Kuwahara  and  Kubota(4] 
extended  the  study  to  low-pressure  combustion  in  order  to  observe  combustion  wave 
structures. 

Kubota[3]  also  investigated  the  combustion  mechanisms  of  nitramine  composite 
propellants  to  determine  controlling  parameters  of  the  burning  rate  in  the  combus¬ 
tion  wave.  He  found  that  the  combustion  wave  structure  was  homogeneous  due  to 


the  melt  layer  formed  by  RDX  and  binder.  He  also  observed  a  preparation  zone, 
similar  to  a  dark  zone  in  the  combustion  of  a  homogeneous  propellant,  whose  height 
strongly  depended  on  the  ambient  pressure.  From  thermocouple  measurements,  he 
concluded  that  the  effect  on  the  burning  rate  of  heat  feedback  from  the  luminous 
flame  was  negligible  in  comparison  to  heat  released  at  the  burning  surface. 

The  effects  of  binders  on  burning  rates  and  thermal  wave  structures  of  HMX 
composite  propellants  were  investigated  by  Kubota[5].  The  thermal  structure  of 
propellants  with  a  high-oxygen  content  binder  was  homogeneous  due  to  the  forma¬ 
tion  of  a  molten  layer  of  HMX/binder  on  the  burning  surface.  The  flame  structure 
of  HMX  propellants  was  found  to  be  a  two-stage  flame  consisting  of  two  reaction 
zones.  The  first-stage  reaction  zone,  which  appears  just  above  the  burning  sur¬ 
face,  is  equivalent  to  the  fizz  zone  of  double-base  propellants;  the  second  stage, 
i.e.,  the  preparation  zone,  is  equivalent  to  the  dark  zone  of  double-base  propellants. 
The  luminous  flame  that  appears  some  distance  above  the  burning  surface  is  pro¬ 
duced  by  the  second  stage  reaction.  The  luminous  flame  and  the  preparation  zone 
were  dis  inguished  by  measurement  on  microphotographs  and  microthermocouples. 
However,  the  thermocouple  outputs  were  not  synchronized  with  the  high-speed  mi¬ 
crophotography. 

The  combustion  of  HMX-based  propellants  was  also  studied  by  Cohen-Nir(7]. 
From  his  cinemicrography  study,  he  argued  that  previous  speculations  regarding 
the  ejection  of  HMX  particles  from  the  burning  surface  were  erroneous. 

Baumann  et.  al[39]  and  Harris[40]  measured  the  temperature  distribution  and 
transient  species  concentrations  in  the  reaction  zone  of  an  RDX-based  composite 
propellant  using  Coherent  Anti-Stokes  Raman  Scattering  (CARS)  to  provide  infor¬ 
mation  for  the  evaluating  the  reaction  mechanism  in  the  decomposition  of  nitramine 
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propellants.  The  species  identified  were  N2, H2,CO, CO2, HCN,  and  NO.  The  pro¬ 
pellant  flame  consisted  of  three  zones:  an  inner  dark  zone,  a  bright  yellow  interme¬ 
diate  zone,  and  an  outer  blue  zone.  The  inner  dark  zone,  which  is  about  4mm  thick 
in  air  at  1  atm,  is  the  reaction  zone  near  the  solid-gas  interface  where  the  surface 
temperature  is  around  900  +  100  K.  In  this  zone,  large  percentages  of  CO  and  H2 
are  formed.  The  final  products  were  identified  as  N2,H2,CO,  and  CO2',  intermedi¬ 
ate  species  such  as  N2O  and  NO2  were  not  observed.  Baumann  and  his  coworkers 
claimed  that  the  intermediate  NO2  at  atmospheric  pressure  was  not  observable 
with  their  CARS  apparatus.  Despite  the  fact  that  no  NO2  species  were  measured, 
they  suggested  that  the  initial  step  of  the  thermal  decomposition  of  nitramines  is 
the  loss  of  NO2  to  form  a  dinitro-RDX  radical,  which  further  undergoes  stepwise 
degradation  to  form  n-nitroformimine  ( H2C  =  N  —  NO 2)  and  HCN.  At  low  tem¬ 
peratures  (below  600  K),  CH2O  and  N2O  form,  whereas  at  higher  temperatures 
(above  700  K),  HCN  and  N2O  form.  Based  upon  the  existence  of  moderate  quanti¬ 
ties  of  HCN,  they  believed  the  high-temperature  pathway  to  be  the  major  route  for 
thermal  decomposition  of  RDX  near  the  surface  of  the  propellant.  This  is  consis¬ 
tent  with  the  decomposition  mechanism  considered  by  Boggs[48]  and  Schroeder[43]. 
They  concluded  that  the  luminous  flame  due  to  conversion  of  NO  to  N2  was  too 
far  away  to  provide  to  the  burning  surface  the  heat  that  determines  the  burning 
rate.  They  believed  that  the  heat  feedback  from  the  gas-phase  reactions  of  RDX 
and  HCN  controlled  the  burning  rate. 

Parr  and  Hanson-Parr[41]  investigated  the  chemistry  and  kinetics  of  the  ig¬ 
nition  and  combustion  of  HMX  using  PLIF (Planar  Laser  Induced  Fluorescence). 
They  measured  species  and  temperature  profiles  during  CO2  laser  ignition  and 
steady-state  deflagration  of  HMX.  They  obtained  time-resolved  PLIF  images  for 
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JVOj,  NO,  CN ,  NH,  H2CO,  OH,  and  OH  rotational  temperatures  during  ig¬ 
nition  and  deflagration  of  HMX.  CN  and  NH  were  found  in  a  delayed  gas-phase 
ignition  kernel  which  rapidly  reformed  into  a  thin  flame  sheet.  They  concluded  that 
NO  and  N02  were  initial  decomposition  products  which  were  consumed  in  a  flame 
in  which  CN  and  NH  were  transient  radical  by-products. 

A  theoretical  combustion  model  of  nitramine  propellants  was  proposed  by  Ben- 
Reuven  et  al.[8-10].  For  the  deflagration  of  RDX,  they  proposed  two  overall  reaction 
steps: 


RDX  -+  -N2  +  N20  +  N02  +  ZCH20  (2  -  3) 

It 

Ich,o  +  no,-,no  +  2-co,  +  3-co^-h,o  (2-4) 

The  proposed  mechanism  included  (1)  the  decomposition  of  RDX  in  the  liquid 
phase,  (2)  the  primary  gas-phase  reaction  following  vaporization  and  gas-phase 
decomposition,  and  (3)  oxidation  of  CH20  by  N02  that  is  the  secondary  gas- 
phase  reaction.  In  the  early  phase  of  their  work  [8],  they  assumed  Eq.(2-2)  to  hold 
for  both  gas-  and  liquid-phase  decomposition  of  RDX.  The  reaction  scheme  was 
modified  later(10]  to  consider  different  reaction  mechanisms  for  the  liquid  phase 
and  gas  phase: 


RDXi  —  3  CH20  +  +  7  NO  +  l  N2 

4  4  8 


(2-4) 


RDX,  —  NO 2  +  N,0  +  3Cff,0  +  ^N, 

It 


(2-5) 


They  found  that  the  near-held  region  with  its  dominant  primary  decomposition  is 
fundamental  in  controlling  heat  feedback  to  the  propellant  surface. 

2.1.3  Scope  of  the  Present  Work 

The  main  objective  of  this  investigation  is  to  develop  a  comprehensive  theory 
for  the  radiative  ignition  of  nitramine  composite  propellants  and  to  evaluate  the 
theory  by  comparison  with  experimental  results  obtained  in  this  investigation.  The 
scope  of  the  present  work  for  theoretical  and  experimental  approaches  is  presented 
below. 


2.1.3. 1  Theoretical  Approach 

The  specific  objectives  of  the  theoretical  study  are: 

1)  to  develop  a  theoretical  model  and  to  obtain  numerical  solutions  for  ignition 
of  nitramine  composite  propellants  under  COa  laser  heating; 

2)  to  validate  the  theoretical  model  with  experimental  data; 

3)  to  study  important  parameters  governing  the  laser  ignition  of  nitramine  com¬ 
posite  propellants. 

2. 1.3.2  Experimental  Approach 

The  specific  objectives  of  the  experimental  study  are: 

1)  to  determine  the  dependence  of  delay  time  for  onset  of  light  emission  (tj,#)  as 
a  function  of  incident  energy  flux  (g"),  ambient  gas  compositions,  and  ambient 
pressure; 

2)  to  observe  the  site  of  initial  light  emission  and  flame  development  processes; 

3)  to  measure  gas-phase  temperatures  in  the  vicinity  of  the  burning  surface; 

4)  to  study  the  surface  structures  of  test  samples  recovered  before  and  after  onset 
of  light  emission;  and 
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5)  to  compare  the  ignition  characteristics  of  a  series  of  nitramine  solid  propellants 
with  those  of  a  conventional  AP/HTPB  propellant. 


2.2  THEORETICAL  APPROACH 


2.2.1  Description  of  Physical  Model 

Physical  processes  involved  in  radiative  ignition  of  composite  solid  propellants 
are  complex.  Figure  1  shows  a  schematic  diagram  of  various  thermophysiochemical 
processes  involved  in  the  radiative  ignition  of  composite  solid  propellants.  In  the 
initial  inert  heating  period,  incident  radiative  energy  heats  up  the  solid  propellant 
material  without  involving  any  chemical  reactions.  Once  the  surface  temperature 
of  the  propellant  reaches  a  certain  point,  a  significant  amount  of  gas  is  generated 
from  the  pyrolysis  of  the  condensed-phase  propellant.  A  melt  layer  may  be  formed 
before  massive  pyrolysis  occurs.  As  the  pyrolyzed  gas  product  is  moving  away  from 
the  propellant  surface,  a  part  of  the  incident  external  radiative  energy  is  absorbed 
by  the  pyrolyzed  gas.  A  portion  of  the  radiative  energy  reaching  the  propellant 
surface  is  reflected  away  from  the  surface.  The  remainder  of  the  radiative  energy 
penetrates  into  the  propellant  and  is  absorbed  by  the  condensed  phase.  Pyrolyzed 
gas  of  the  RDX  filler  may  undergo  a  primary  monoflame.  The  gas  product  of 
the  primary  flame  diffuses  and  reacts  with  the  diffusing  decomposed  gas  from  the 
binder.  Diffusion  of  external  oxidizing  gas  may  accelerate  the  gas-phase  chemical 
reactions. 

The  heterogeneous  structure  of  nitramine  composite  propellants  is  shown  in 
Fig.  2.  The  physical  model  of  the  present  study  considers  an  RDX  filler  particle 
embedded  in  the  binder  matrix.  The  shape  of  the  filler  particles  is  approximated  by 
a  cylindrical  pellet  of  radius  R\  and  thickness  L,  enveloped  by  a  binder  of  radius  R.2- 
Figure  3  shows  the  statistically  averaged  element  and  finite  difference  grid  pattern. 
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Fig.  1.  Physical  and  Chemical  Processes  in  Composite  Solid  Propellant 
Ignition  under  Radiative  Heating 
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2,2.2  Assumptions 


The  following  assumptions  are  employed  in  the  mathematical  model: 

1)  The  solid  propellant  and  gas  phase  are  two-dimensional  and  axisymmetric. 

2)  The  gas  pressure  is  constant. 

3)  The  gas-phase  chemical  reactions  can  be  described  by  second-order  Arrhenius 
kinetics. 

4)  The  pyrolysis  processes  are  expressed  by  single  zero-order  reactions  of  the  Ar¬ 
rhenius  Law. 

5)  The  gas  mixture  obeys  the  perfect  gas  law. 

6)  The  specific  heat  of  all  gaseous  species  is  assumed  to  be  the  same  and  constant. 

7)  The  binary  diffusion  coefficients  of  all  gaseous  species  are  equal. 

8)  The  incident  radiation  is  absorbed  in-depth  by  the  solid  phase  satisfying  Beer’s 
Law. 

9)  The  radial  bulk  velocity  is  much  smaller  than  the  axial. 

10)  The  thermal  properties  of  the  solid  are  independent  of  temperature. 

11)  The  Lewis  number  is  unity. 

The  assumption  of  a  two-dimensional  geometry  closely  represents  a  composite 
propellant,  while  being  mathematically  tractable.  Since  there  is  no  external  driving 
force  for  the  ambient  gas  for  radiative  ignition,  the  pressure  of  the  ambient  gas  is 
nearly  constant.  This  assumption  eliminates  the  momentum  equations.  The  as¬ 
sumptions  for  the  gas-phase  chemical  reactions  are  due  to  the  lack  of  knowledge 
of  the  complex  chemical  reactions  in  the  gas  phase  and  the  prohibitive  computa¬ 
tional  time  involved.  Further  discussion  will  be  presented  in  subsequent  sections, 
and  as  further  information  is  gathered,  the  more  refined  reaction  schemes  will  be 
incorporated.  The  pyrolysis  mechanisms  of  RDX  filler  remain  controversial  even 


after  extensive  studies  of  the  mechanisms.  This  subject  will  be  discussed  further  in 


a  later  section.  Assumptions  (5),  (7),  and  (11)  are  standard  in  combustion  studies. 
Assumptions  (6)  and  (10)  are  made  to  simplify  the  problem.  Assumption  (9)  holds 


true  for  radiative  ignition. 


The  mathematical  model  consists  of  governing  equations  for  the  solid-phase 


and  gas-phase  regions.  In  the  solid  phase,  two  energy  equations  are  considered; 


one  for  the  RDX  filler  and  the  other  for  the  binder.  Source  terms  in  these  equa¬ 


tions  include  contributions  due  to  in-depth  radiation  absorption  and  pyrolysis.  The 


energy  equations  are  coupled  to  the  gas-phase  conservation  equations  through  the 
heat-flux  balance  at  the  solid-gas  interface.  Gas-phase  behavior  is  described  by 


mass,  energy,  and  species  conservation  equations.  For  the  RDX-based  composite 


propellants,  five  different  species  ( CH?Q ,  NO2,  RDX  vapor,  gaseous  binder,  and 
products)  are  considered  to  be  present  in  the  gas  phase. 


RDX  filler: 


,3T  rhR  3T.  d  (,  dT\  1  d  /.  dT\ 
m  j  dz  )  -  gt  (  *•'  3*  )  +  r  Sr  (*K'*r  dr)+  +  ,R'! 


(2-6) 


Binder: 


Pfl1#CBl«(  dt  +  dz)  dz  B,a  dz  )  +  rdr{  B,a  ^  J  +  ?B,rad  +  <1b>p  yrc 

(2  -  7) 


where 


Cad  =  -^=^/o(l-Pr,,V^ 


ft; 

■ft; 

ft 


Vi 

J»Sh*L. 
•  !» 

It*. V 


Lv;: 


Avi  V.SsV«*v 


iilpyro  =  ~Zj,PVroQ},pvroCZp(-EjlPyro/ RhT), 


and 


j  =  R{RDX)  or  B  (binder)  (2-8) 


Gas- Phase  Equations 


The  gas-phase  conservation  equations  are 
Continuity: 


dpg  d  .  . 

~m  +  ai  ^  =  0 


(2-9) 


Energy: 


„  dT  „  dT  d  f.  dT\  1  d  /.  dT\  tn 

*c’dt  +  -  ai  (*»  57 j  +  ;  *  (A'jrd7 )  +  9«  (2 “ 10) 


Species: 


dYj  dY,  d  (  ndYj\  1  d  (  .  /rt 

+ =  ai  (,,'p  ar)  +  (','Dra^J  +w<  (2 "  n) 


where  j  =  1,  2,  3  or  4  for  the  gaseous  species  C/fjO,  N^O,  RDX  vapor,  and 
binder  pyrolysis  gas,  respectively. 

The  source  terms  q'g"  and  a/'"  are  obtained  by  using  chemical  kinetics  informa¬ 
tion  described  in  a  later  section.  The  equation  of  state  for  the  gas  phase  is 


PMW 
Pt  “  RUT 


(2  -  12) 
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Initial  and  Boundary  Conditiona 

The  initial  condition  for  the  solid  phase  (both  RDX  and  binder)  is 

at  t-  0  :  r(0,r,*)  =  T,  (2-13) 

Boundary  conditions  for  the  solid  phase  are  as  follows.  The  temperature  far 
from  the  interface  is  the  same  as  the  propellant  initial  temperature,  i.e., 


at  m  -*  -oo  :  T  =  Ti 


(2  -  14) 


Temperature  continuity  at  various  interface  requires 


rL*f-  rU*r 


and 


=  t  |„_t. 

On  symmetric  surfaces,  adiabatic  conditions  are 


(2  -  15) 


(2  -  16) 


dT 

on  r  =  0  :  —  =  0 

or 


(2  -  17) 


dT 

on  r  =  Ri  :  —  =  0 
dr 

At  solid-gas  and  RDX-binder  interfaces,  the  heat  flux  balance  gives 


(2  -  18) 


dT  dT 

at  z  =  Zt-g  :  _  =  kg~dg\*,  „+  +  cts-gl  +  rt ,PtT(cPi3  —  cPtg)  +  qa_g 

(2  -  19) 

where  g"  is  the  net  heat  generation  at  the  interface. 
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at  m  =  —L  = 


-t.  +  ai.-ar./.e-'’*1'  (2-20) 


a»r 

«t  r  =  and  -£<*<0  :  (2  ~  21) 

Initial  conditions  for  the  gas-phase  equation  (for  j  =  1,  2,  3,  and  4)  are 


at  t  =  0  :  v„(0,  r,  m)  =  0;  r(0,  r,z)  =  r#1»;andy(0,r,s)  =  Yiti  (2-22) 

Boundary  conditions  for  the  gas  phase  are  as  follows.  The  symmetric  conditions 
at  the  centerline  and  outer  boundary  give 


dT  n  dYj 

at  r  =  0  ~  =0  ; 


(2  -  23) 


„  dT  n  dYj 
on  r  =  ft  ■.  —  = 


Far  away  from  the  surface, 


(2-24) 


dT  „  dYj  „ 
*  -  00  :  aT=0i  ~dt*° 

The  temperature  continuity  at  the  interface  gives 


(2  -  25) 


on  e  =  z,-g  :  rl  =  T 


(2  -  26) 


The  overall  and  individual  species  mass-flux  balance  at  the  solid-gas  interface 
is  given  as  follows. 


ivsvartiKvaviva  vtraiMiUt  vRSftkrtimvvu-i  tra  JV  \r%  L**  vis  ■.'('.’v  L-**  k.*n  W>  \Tm 


»  *  K  4 

V*. 


on  z  =  z,-g  :  pgv„  =  parbt 


(2  -  27) 


and 


p°v‘Yi\'.-,- =  ~  p*v~d  ~  (2 " 28) 

2. 2. 3. 2  Phvsiochemical  Consideration  for  Pvrolvsis  and  Ignition 

Description  of  Physical  and  Chemical  Processes 

Physical  and  chemical  processes  involved  in  the  ignition  and  combustion  of 
nitramine  composite  propellant  are  described  below.  Figure  4  shows  a  general  dia¬ 
gram  compiled  from  various  studies  of  nitramine  propellant  decomposition,  ignition, 
and  combustion.  It  is  evident  that  the  ignition  and  combustion  of  nitramine  com¬ 
posite  propellants  involves  many  complex  physical  and  chemical  processes.  When 
heated  by  ambient  gases,  the  solid-phase  nitramine  crystal  and  binder  undergo  inert 
heating  in  the  subsurface  region.  Near  the  surface  of  the  nitramine  crystal,  phase 
transitions  may  occur.  Decomposed  radicals  could  be  trapped  between  molecules  or 
within  molecules.  This  produces  the  so-called  “Cage  Effect”  proposed  by  Fifer[49]. 
Trapped  radicals  can  recombine  due  to  their  close  proximity  and  low  mobility.  This 
recombination  process  is  generally  believed  to  be  the  reason  for  the  higher  acti¬ 
vation  energy  of  solid-phase  decomposition  than  that  of  liquid-phase  or  gas-phase 
decompositions.  It  was  reported  by  Brill[50]  that  the  phase  transition  has  an  acti¬ 
vation  energy  similar  to  that  of  solid-phase  decomposition.  He  postulated  that  the 
rate  ontrolling  step  of  solid-phase  thermal  decomposition  is  due  to  disruption  of 
intei  olecular  forces. 
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At  the  interface  between  the  gas  phase  and  condensed  phase,  a  foam  or  liquid  melt 
layer  could  form  above  the  solid-phase  crystal  and  binder.  The  foam  layer  may 
contain  various  sizes  of  bubbles.  The  melt  layer  of  the  nitramine  filler  could  form  at 
a  temperature  lower  than  the  melting  point  of  pure  nitramine.  This  phenomenon 
was  called  "liquefaction”  by  Boggs [48].  Boggs  stressed  that  liquefaction  is  not  a  true 
melting  process.  A  melt  of  nitramine  filler  and  binder  could  be  mixed  to  fo’-m  an 
energetic  melt  layer.  Vaporization  of  the  liquid  melt  can  occur  on  the  surface  of  the 
melt  layer.  In  the  meantime,  thermal  decomposition  of  the  nitramine  filler  in  the 
liquid  phase  can  yield  decomposed  radicals  and  fragmented  nitramine  molecules  in 
the  initial  decomposition  step.  Further  decomposition  of  the  gaseous  species  formed 
at  the  initial  step  and  vapor-phase  decomposition  will  yield  slightly  oxidizer-rich 
gases.  Some  of  the  gaseous  species  from  the  oxidizer-rich  gases  could  attack  the 
surface  of  the  melt  layer.  Chemical  species  of  the  oxidizer-rich  gases  will  react  with 
each  other  exothermically  to  form  the  first-stage  reaction  zone.  Before  reacting  in 
the  final  flame,  the  gaseous  products  from  the  first-stage  reaction  zone  may  pass 
a  preparation  zone  where  reaction  kinetics  are  retarded.  The  gaseous  species  from 
the  preparation  zone  and  the  decomposed  fuel  gases  will  react  eventually  to  form 
the  final  luminous  flame. 

Radiative  heat  could  be  transferred  from  the  luminous  flame  to  the  first  reaction 
zone  and  the  foam  or  melt  layer.  Simultaneously,  conductive  heat  feedback  to  the 
condensed-phase  surface  from  the  first  reaction  zone  will  occur.  The  mechanisms 
depicted  in  Fig.  4  are  based  upon  observations  and  a  physical  understanding  of  the 
ignition  and  combustion  processes. 

Thermal  Decomposition  of  Nitramine 

Since  it  has  been  considered  that  the  overall  chemical  processes  in  the  deftagra- 


28 


tion  of  nitramine  composite  propellants  are  dominated  by  thermal  decomposition  of 
pure  nitramine,  extensive  studies  have  been  conducted  on  the  subject.  Schroeder[43- 
47],  Boggs[48],  Fifer[49],  and  Dubovitsii  and  Korsunskii[58j  reviewed  the  literature 
on  the  kinetics  and  reaction  mechanisms  of  the  thermal  decomposition  of  nitramines, 
especially  RDX  and  HMX.  Despite  extensive  studies  on  the  thermal  decomposition 
of  RDX  and  HMX,  the  detailed  mechanisms  are  generally  unknown.  Even  the  initial 
step(s)  of  the  reaction  is  controversial. 

Various  mechanisms  proposed  have  been  critically  reviewed  by  Schroeder[43- 
47].  An  excellent  arrangement  of  the  probable  mechanisms  has  been  prepared  by 
Shaw  and  Walker [59].  The  main  question  concerning  the  primary  step  (or  steps)  of 
thermal  decomposition  of  RDX  and  HMX  is  whether  it  is  N  —  N  bond  fission  or 
C  —  N  bond  cleavage.  Fission  of  N  —  N02  bond  has  generally  been  considered  as  an 
initial  step  of  the  decomposition  mechanism  of  RDX  and  HMX.  Growing  evidence 
has  shown  that  C  —  N  bond  cleavage  occurs  first,  identified  mostly  by  mass  spec¬ 
troscopic  (>50, 61]  and  labelling  studies  accompanying  rapid  pyrolysis[62].  The  initial 
decomposition  schemes  are  (1)  elimination  of  CH2NNO2,  (2)  N  —  N02  fission,  (3) 
homolytic  C  —  N  cleavage,  (4)  HONO  elimination,  (5)  concerted  depolymerization 
to  4  CH2NNO2 ,  and  (6)  transfer  of  an  oxygen  atom  from  an  —N02  group  to  a 
neighboring  group. 

A  mechanism  proposed  by  Karpowicz,  Gelfand,  and  Brill[50]  is  quite  different 
from  those  of  covalent  bond  cleavage.  Based  on  observations  of  their  phase  tran¬ 
sition  experiments,  it  was  proposed  that  disruption  of  intermoiecular  electrostatic 
forces  (rather  than  covalent  bond  cleavage)  controls  the  rate  of  thermal  decomposi¬ 
tion  of  HMX.  As  mentioned  above,  Fifer  also  postulated  the  "Cage  Effect,"  which 
represents  an  induced  recombination  of  radicals  generated  in  the  initial  dissociation 


step.  He  suggested  that  cage  escape  (the  escape  of  decomposed  gaseous  species 
from  the  condensed  region)  might  be  rate-controlling  for  the  condensed-phase  de¬ 
composition.  In  general,  both  Fifer  and  Brill  consider  intermolecular  actions  to  be 
dominant  factors  in  condensed-phase  decomposition. 

Kinetic  Model 

Based  on  the  information  discussed  above,  a  kinetic  model  was  constructed. 
The  reaction  mechanisms  used  in  the  present  model  can  be  extended  to  include 
various  initial  decomposition  mechanisms  discussed  previously.  The  initial  step  can 
be  N  —  NO2  fission,  C  —  N  bond  cleavage,  HONO  elimination,  or  other  mecha¬ 
nisms.  The  present  investigation  considers  the  following  five  representative  chemical 
reactions:  (1)  vaporization  of  nitramine  filler  followed  by  thermal  decomposition  to 
form  oxidizer-rich  gases,  (2)  initial  thermal  decomposition  of  liquid  phase  nitramine 
to  form  decomposed  radicals  ( DR )  and  fragmented  RDX  ( RDX ),  (3)  further  de¬ 
composition  of  RDX •  to  form  intermediate  oxidizer-rich  gases  which  then  react  in 
the  first-stage  reaction  to  form  Oxg,  (4)  binder  pyrolysis  to  form  fuel-rich  gases, 
and  (5)  reaction  of  Oxg  and  Fg  to  form  final  products  in  luminous  flame: 


RDXi - - *  RDXg- - t - ♦  Oxg 

vaporization  decomposition 


(2  -  29) 


RDXi  thermal  decomposition  ^Rg  +  RDXg 


(2  -  30) 


DR'g  +  RDXg 


^3o 

further  decomposition 


Intermediate  Oxidizerrich  Gases 


k^b 


first  stage  reaction 


Ola 


(2  -  31) 
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Oxg  +  Fg  — r 


k$ 


-*■  products g 


(2  -  33) 


9  9  fined  flame 

Depending  upon  the  options  considered  for  radical  initiation,  DR •  could  repre¬ 
sent  NO’2  ,C  H2NNO’2  ,  HO  NO'  or  any  other  reactive  species.  Some  of  the  kinetic 
data  required  as  input  for  the  thermal  decomposition  of  nitramine  filler  is  available 
from  various  decomposition  experiments.  Other  kinetic  constants  such  as  Arc  re¬ 
quired  for  the  above  reaction  model  must  be  based  upon  the  kinetic  data  of  similar 
reactions. 


2.2.3.3  Modeling  of  Source  Terms 

In  order  to  solve  the  governing  equations,  source  terms  (q",  qg",  w",  and  u>"') 
must  be  expressed  as  functions  of  other  variables  and  known  parameters.  In  model¬ 
ing  the  source  terms,  it  is  implicitly  assumed  that  the  rate-controlling  reaction  rate 
constants  have  Arrhenius  dependence  on  local  temperature  and  are  given  by 

k  =  zTme~ElR'T  (2-34) 

The  gas-phase  mass  fractions,  Yi,  Yj,  Y3,  Y4,  and  Y$,  represent  the  RDX  vapor, 
NOi,  CH20 ,  binder  decomposed  gas,  and  products,  respectively.  The  source  terms 
can  be  expressed  as  follows: 

<  =  Z3 ^-pgTexp(-E3/RvT)  -  Z6^-pgYiY4cxp(-E5/ RUT)  (2  -  35) 


=  Z3±-p2gY2Y3Texp(-E3/RuT) 


(2  -  36) 
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=  K'Affa  +  <4«l)  (2  -  39) 

where  Alfs  is  the  heat  of  reaction  per  unit  mass  of  RDX  for  the  reaction  of  Eq.(2- 
31),  and  AJSs  is  the  heat  of  reaction  per  unit  mass  of  fuel  for  the  reaction  of 
Eq.(2-33). 


2.2.3.4  Transformation  of  Governing  Equations 

In  view  of  the  large  temperature  and  species  concentration  gradients  near  the 
propellant  surface,  a  coordinate  transformation  is  used  to  obtain  an  exponential 
finer  grid  size  in  the  transformed  coordinates.  The  following  exponential  transfor¬ 
mation  is  used  for  the  axial  coordinate. 

Solid  phase: 

s  =  exp(AaZ)  -  1  or  Z  =  +  s),  (-00  <  Z  <  0)  (2  —  40) 

Gas  phase: 

a  =  1  -  exp(-AgZ)  or  Z  =  —  «)» (0  <  Z  <  oo)  (2  —  41) 

"  Ag 

This  transforms  £  =  -op,  0,  oo  to  s  --1,  0,  l,  respectiveiy.-IhgradigiToordinates 
remain  unchanged. 
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2.2.4  Numerical  Method 


An  implicit  finite  difference  scheme  was  used  to  solve  the  governing  equations. 
A  central  difference  scheme  was  used  to  approximate  axial  derivatives.  A  three- 
point,  variable  mesh,  Allen’s  method  was  employed  to  represent  radial  derivatives, 
because  the  conventional  central-difference  scheme  does  not  accurately  represent 
radial  derivatives  when  the  radial  coordinate  is  very  small.  Time  derivatives  were 
approximated  by  a  generalized  Crank-Nicolson  scheme.  A  quasilinearization  tech¬ 
nique  was  used  to  linearize  the  inhomogeneous  terms  of  governing  equations.  The 
resulting  set  of  simultaneous  finite-difference  equations  are  being  solved  by  a  suc¬ 
cessive  overrelaxation  iterative  scheme. 


2.3  EXPERIMENTAL  APPROACH 


2.3.1  Introduction 

An  experimental  investigation  of  the  laser  ignition  of  nitramine  composite  pro¬ 
pellants  was  initiated  to  explore  the  physical  and  chemical  processes  involved  in 
interactions  between  an  incident  laser  beam  and  composite  propellants.  Observa¬ 
tions  and  understandings  gained  from  these  experiments  are  useful  in  developing 
new  igniter  systems,  and  are  very  important  for  validating  proposed  ignition  models. 
Thus,  the  experimental  approach  was  carefully  selected,  baaed  on  the  possibility  of 
the  comparison  with  theoretical  results.  Knowledge  obtained  from  this  radiative  ig¬ 
nition  can  be  extended  to  the  ignition  processes  under  other  types  of  energy  transfer 
modes,  either  convective,  conductive,  or  combinations  of  the  three  different  types 
of  energy  transfer  modes. 

In  this  chapter,  the  experimental  apparatus  and  their  operations  will  be  dis¬ 
cussed  in  the  first  section.  Then,  the  data  acquisition  system  will  be  discussed, 
followed  by  the  description  of  the  propellant  samples  used  in  the  experiment.  The 
laser-beam  intensity  distribution  which  is  important  in  defining  the  energy  flux  to 
the  sample  surface  is  discussed  next.  The  examination  of  beam-intensity  distri¬ 
bution  will  be  followed  by  description  of  the  techniques  used  in  measuring  both 
temperature  and  gaseous  species.  Finally,  the  schlieren  technique  used  for  visual¬ 
ization  of  gasification  and  flame  development  will  be  discussed. 

2.3.2  Experimental  Apparatus  and  Operation 

The  schematic  diagram  of  the  CO2  laser  ignition  test  setup  is  shown  in  Fig.  5. 
A  low-pressure  ignition  test  chamber,  250  x  250  mm  and  250  mm  high,  is  mounted 
on  an  X-Y  positioning  bed  to  provide  convenience  in  aligning  the  laser  beam  on 
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C02  LASER 
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the  teat  sample.  The  propellant  sample  is  mounted  on  a  black-anodised  aluminum 
sample  holder  located  . the  center  of  the  test  chamber.  In  order  to  minimise 
pressure  change  during  the  ignition  test,  the  sise  of  the  test  chamber  was  designed 
to  be  large  enough  so  that  the  effect  of  the  product  gas  generated  from  the  sample 
on  the  chamber  pressure  during  the  test  would  be  negligible.  The  gas  supply  line  to 
provide  different  gases,  and  the  purge  line  to  remove  the  combustion  gas  products, 
are  connected  to  the  stainless  steel  top  cover  of  the  test  chamber. 

The  four  side  walls  of  the  low-pressure  chamber,  suitable  for  a  pressure  range 
of  10~3  to  2  atm,  were  fabricated  with  plexiglass  (Fig.  6).  The  transparent  side 
walls  of  the  chamber  provide  full  access  to  ignition  events  from  all  directions. 
The  feedthroughs  for  the  photodiode,  thermocouples,  and  gas-sampling  probes  are 
mounted  on  the  side  walls  of  the  chamber.  For  schlieren  photography,  a  pair  of 
schlieren  quality  windows  were  mounted  on  the  two  opposite  side  walls  of  the  test 
chamber.  The  C02  laser  beam  reaches  the  propellant  sample  by  passing  through 
the  KC1  window  mounted  on  top  of  the  chamber.  A  gas  control  system  (Fig.  7) 
was  also  designed  and  constructed  to  provide  the  desired  initial  gas  composition 
and  chamber  pressure. 

A  high-powered  Coherent  Super  48  C02  laser  (Fig.  8)  h  used  as  a  radiant 
energy  source  for  the  ignition  study  of  nitramine  composite  propellants.  The  maxi¬ 
mum  power  of  the  C02  laser  is  800  W  in  the  CW  mode,  and  3500  W  in  the  pulsed 
mode.  The  laser  time  can  be  controlled  precisely  by  preselecting  the  time  period 
of  irradiation.  The  C02  laser  generates  a  signal  to  trigger  all  the  data  acquisition 
systems. 

The  C02  laser  beam  generated  from  the  laser  system  travels  through  a  beam 
delivery  system  consisting  of  a  series  of  silicon  mirrors.  The  laser  beam  is  confined 
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b)  Close-up  View 
Fig.  6.  Low  Pressure  Chamber 
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Fig.  7  A  Schematic  Diagram  of  Gas  Control  System 
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within  the  anodized  beam  guiding  tubings.  Since  the  CO2  laser  beam  is  invisible 
(wave  length  =  10.6  /tm),  a  beam  combiner  was  designed  and  fabricated  to  combine 
the  invisible  CO2  laser  beam  with  an  He-Ne  laser  beam  so  that  the  two  beams  travel 
together.  One  of  the  major  advantages  in  combining  an  He-Ne  laser  with  a  CO2 
laser  is  to  provide  convenience  in  beam  alignment. 

2.3.3  Data  Acquisition  and  Analysis  Systems 

Three  different  types  of  cameras  (a  high-speed  video  camera,  a  high-speed 
movie  camera,  and  a  conventional  video  camera)  are  used  to  observe  flame  initiation 
and  development.  A  Spin-Physics  2000  high-speed  video  camera  (12,000  pps  max.) 
is  used  to  film  flame  development.  A  quick  review  of  the  recorded  ignition  events  is 
possible  with  this  high-speed  video  system,  and  the  recorded  images  can  be  analyzed 
on  the  monitor  by  moving  reticle  lines.  Data  can  be  transferred  to  an  advanced 
Qauntex  image  analyzer  for  processing  and  storage.  A  high-speed  movie  camera 
(HYCAM  44,000  pps  max.)  is  used  to  obtain  a  more  detailed  record  of  gas  evolution 
and  flame  development  during  pyrolysis  and  ignition.  A  conventional  video  camera 
rystem,  despite  its  low  filming  rate  (ca.  25  pps),  is  also  useful  in  observing  the 
overall  event  in  relatively  good  spatial  resolution.  In  addition  to  the  movie  cameras, 
a  silicon  photodiode  is  used  to  detect  the  onset  of  first  light  emission. 

Data  from  the  thermocouples  and  the  photodiode  are  recorded  on  a  Nicolet 
digital  oscilloscope.  The  data  recorded  by  the  Nicolet  oscilloscope  is  stored  on  floppy 
diskettes,  and  then  transferred  to  an  IBM- AT  microcomputer.  The  transferred  data 
is  processed  and  analyzed  by  the  microcomputer  and  plotted  on  a  laser-jet  printer. 


The  composition  and  thermophysical  properties  of  the  ingredients  of  the  pro¬ 
pellant  samples  used  in  this  study  are  given  in  Tables  I  and  II.  Three  different 
types  of  RDX-based  LOVA  propellants  (BLX4,  BLX8,  and  BLX9)  processed  at 
the  Naval  Weapons  Center(NWO)  are  tested.  (BLX  stands  for  baseline  explosive.) 
Both  BLX4  and  BLX9  propellants  contain  RDX  fillers  and  GAP  binder,  but  with 
different  energetic  plasticizers  (TMETN  and  BTTN,  respectively).  BLX8  propel¬ 
lant  contains  RDX  filler  and  GAP  binder,  but  no  energetic  plasticizers.  The  RDX 
particles  in  the  BLX-series  propellants  are  divided  in  two  different  classes:  classes  A 
and  E.  The  average  size  of  the  particle  in  the  class  E  is  4/*m,  while  class  A  contains 
a  wide  range  of  particle-size  distribution.  In  order  to  compare  the  ignition  behav¬ 
ior  of  BLX-series  of  nitramine  propellants  with  conventional  AP-based  composite 
propellants,  a  propellant  processed  at  Atlantic  Research  Corporation  (ARC)  with 
73%  AP  (average  particle  size  —  20  urn)  and  27%  HTPB  was  also  used. 


In  order  to  define  the  incident  energy  flux  to  the  sample  surface,  a  uniform 
distribution  of  laser-beam  intensity  is  desired.  The  intensity  distribution  of  the 
COj  laser  beam  was  examined  by  using  a  fine  wire  thermocouple  and  a  pyroelectric 
detector.  A  fine  wire  thermocouple  was  mounted  on  a  pole  attached  to  an  x-y  po¬ 
sitioning  bed.  The  thermocouple  travelled  through  the  laser  beam  while  measuring 
the  intensities  of  the  beam  at  different  equally  spaced  locations.  Since  the  exposure 
time  of  the  thermocouple  to  the  beam  was  sufficiently  short,  the  slopes  of  the  ther¬ 
mocouple  signals  were  quite  straight,  indicating  that  the  convective  cooling  effect 
of  the  thermocouple  was  negligible.  The  relative  magnitude  of  the  incident  energy 
flux  was  determined  based  on  the  slopes  of  the  signals.  Since  the  heat-transfer 


Table  I.  Composition  of  Test  Propellants 
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mechanism  from  the  laser  to  the  thermocouple  is  quite  complicated,  no  attempt 
was  made  to  determine  the  absolute  amount  of  the  energy  flux.  A  pyroelectric 
detector(PZT  350),  in  conjunction  with  a  chopper  whose  frequency  was  100  Hz,  was 
used  to  verify  the  intensity  distribution  measured  by  the  fine-wire  thermocouple. 
The  profile  of  beam  intensity  distribution  showed  that  a  low  intensity  region  exists 
near  the  center  of  the  beam. 

Two  different  approaches  were  taken  to  achieve  uniformity  of  beam  intensity 
distribution.  The  combination  of  a  Kaleidoscope  and  a  rotating  mirror  was  one 
method  used  to  average  the  original  laser  beam.  The  other  technique  was  to  modify 
the  intensity  distribution  with  a  mask  made  of  thin  stainless  steel.  The  latter 
method  proved  to  be  more  convenient  and  effective. 

The  intensity  level  of  the  laser  beam  was  also  examined  using  the  pyroelectric 
detector.  An  initial  transient  peak  was  observed.  However,  the  delay  time  and  the 
magnitude  of  the  initial  transient  was  not  considered  significant  to  the  ignition  test 
results. 

2.3.6  Temperature  Measurement 

The  surface  temperature  and  thermal  wave  structure  in  the  gas  phase  during 
the  ignition  of  propellants  are  very  important  in  understanding  the  mechanism  of 
pyrolysis  and  ignition.  A  fine  wire  thermocouple  technique  was  used  to  measure 
the  temperatures  of  the  propellant  surface  and  the  gas  evolved  from  the  sample 
surface  during  laser  irradiation.  Figure  9  shows  the  configuration  of  the  thermo¬ 
couples  mounted  with  a  propellant  sample.  The  fine  wire  thermocouple  (dia.  25  /xm 
Pt-Pt/13%Rh)  was  supported  by  thicker  extension  wires.  Since  the  gas  products 
escaping  from  the  sample  surface  can  push  the  thermocouples  upward,  lead  weights 
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Fig.  9  Configuration  of  Fine  Wire  Thermocouple  Mounting 
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Fig.  11  Schematic  Diagram  of  Temperature  Measurement 


were  suspended  on  the  thermocouple  wires  to  maintain  the  location  of  the  thermo¬ 
couple  junctions.  Figure  10  shows  the  magnified  view  of  the  thermocouple  junction. 
In  order  to  minimize  the  size  of  thermocouple  bead,  a  constant  current  technique 
was  used  to  weld  the  thermocouple  wires.  Signals  generated  from  the  thermocou¬ 
ples  were  recorded  on  a  Nicoiet  oscilloscope  through  an  amplifier.  Data  recorded 
by  the  scope  was  processed  and  analyzed  by  an  IBM-AT  microcomputer.  Figure  11 
shows  the  schematic  diagram  of  the  temperature  recording  and  processing  system. 

2.3.7  Gas  Species  Measurement 

Analysis  of  the  composition  of  the  gas  products  from  the  pyrolysis  and  igni¬ 
tion  of  solid  propellants  is  important  to  understanding  the  reaction  mechanisms 
governing  the  processes.  A  Varian  3700  Gas  Chromatograph  was  used  to  analyze 
the  gaseous  species  generated  from  the  sample  during  pyrolysis  and  ignition.  Gas 
samples  were  taken  through  the  gas  sampling  port.  The  stable  species  from  the 
product  gases  were  detected  by  using  two  different  types  of  columns:  a  molecular 
sieve  (MS  5A)  and  a  series  of  Porapak  Q  and  R.  In  order  to  minimize  further  reac¬ 
tions  of  the  initial  decomposed  gases,  low  pressure  conditions  were  chosen  for  the 
pyrolysis  study.  A  GC/MS  system  was  also  used  to  :onfirm  the  peaks  observed  in 
the  gas  chromatography. 

2.3.8  Surface  Structure  Observation 

Scanning  electron  micrography  was  used  to  examine  the  surface  structures  of 
propellant  samples  before  and  after  exposure  to  COj  laser  radiation.  The  microg¬ 
raphy  provided  important  clues  in  determining  the  relative  response  of  the  RDX 
filler  and  the  binder  systems  to  the  incident  laser  beam.  It  was  also  important  in 
ascertaining  the  existence  of  an  energetic  melt  layer  during  ignition  processes. 
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2.3.9  Schiieren  Technique 


One  of  the  common  visualisation  techniques  for  the  study  of  combustion  phe¬ 
nomena  is  schiieren  photography.  Observations  of  the  detailed  gas  evolution  and 
mixing  are  possible  using  a  schiieren  optical  system  during  pyrolysis  and  ignition  of 
the  propellant  samples.  Figure  12  shows  the  configuration  of  the  schiieren  system 
used  in  this  study.  The  "Z"  type  schiieren  optical  system  employs  two  parabolic 
first-surface  mirrors,  each  with  a  diameter  of  6  inches.  The  light  source  is  a  con- 
tinous  tungsten-halogen  lamp.  When  a  Spin-Physics  high-speed  video  camera  is 
used  to  record  the  events,  a  monochromatic  visualisation  technique  is  used.  A  color 
schiieren  technique  developed  by  Professor  Gary  S.  Settles  was  used  to  take  pictures 
with  a  high  speed  movie  camera. 


Fig.  12  Schlieren  Optical  System 


2.4  RESULTS  AND  DISCUSSION 


2.4.1  Experimental  Results 

In  this  section,  the  delay  times  for  onset  of  light  emission  are  presented  and 
discussed.  The  effect  of  oxygen  concentration  in  the  ambient  gas  on  delay  times  is 
also  discussed.  Gas  evolution,  flame  initiation,  and  flame  development  are  examined 
from  direct  and  the  schlieren  pictures  recorded  by  high-speed  video  camera,  high¬ 
speed  movie  camera,  and  conventional  video  camera.  Surface  temperatures  and  the 
thermal  wave  structure  of  the  gas  phase  are  presented  and  discussed  to  explain  the 
thermal  processes  that  occur  during  pyrolysis  and  ignition  of  nitramine  composite 
propellants.  Analysis  of  product  gases  generated  from  the  test  propellants  can 
reveal  the  mechanisms  of  chemical  reactions  involved  in  the  pyrolysis  and  ignition. 
Results  of  the  gas  analysis  by  a  gas  chromatography  are  discussed.  Examination 
of  the  optical  and  scanning  electron  microscope(SEM)  pictures  of  the  propellant 
surfaces  before  and  after  test  illustrates  the  effects  of  incident  energy  fluxes. 

2.4,1. i  Ignition  Delay  Time 

Delay  times  for  the  onset  of  light  emission  in  the  gas  phase  were  determined 
from  photodiode  signals.  The  data  obtained  by  the  photodiode  were  compared 
with  the  delay  times  estimated  from  analysis  of  the  video  pictures  recorded  by  the 
Spin-Physics  high-speed  video  camera. 

Figure  13  shows  the  measured  delay  times  of  the  BLX4  propellant  as  a  func¬ 
tion  of  incident  laser  energy  fluxes  for  three  different  oxygen  concentrations  in  the 
ambient  gas.  The  data  for  1%,  596,  and  2196  O2  show  the  same  trends:  the  delay 
time  decreases  as  the  energy  flux  increases.  For  21$6  O2,  the  delay  times  are  shorter 
than  those  at  lower  oxygen  concentrations.  This  indicates  that  the  oxygen  in  the 
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ambient  gas  participates  in  the  chemical  reactions  that  produce  a  luminous  flame. 
The  slope  of  the  curve  for  the  21%  O2  case  is  less  steep  than  those  for  the  1%  and 
5%  cases.  As  the  incident  energy  increases,  the  effect  of  the  oxygen  concentration  on 
delay  time  becomes  less  pronounced.  For  the  lower  oxygen  concentration  case(l% 
and  5%  Oa),  no  luminous  flame  in  the  gas  phase  was  observed  below  about  600 
VF/ema. 

The  delay  times  for  BLX8  are  given  in  Fig.  14.  For  5%  O3,  a  luminous 
flame  was  observed  in  a  narrow  range  of  the  energy  flux.  At  the  lower  oxygen 
concentration  (1%  O 2),  a  luminous  flame  was  barely  visible.  For  21%  O 2,  a  slope 
break  was  found  at  the  incident  energy  flux  of  about  700  W/cm2.  It  is  believed 
that  the  mechanisms  dominating  the  gas-phase  reactions  are  changing  across  that 
boundary.  The  effect  of  the  absorption  of  the  incident  laser  energy  by  the  gaseous 
products  may  become  significant  at  higher  energy  fluxes. 

Figure  15  shows  the  delay  times  of  the  BLX9  propellant  samples.  The  delay 
times  of  BLX9  are  quite  similar  to  those  of  BLX4.  The  slope  of  the  delay  times  of 
BLX9  at  21%  O3  is  slightly  flatter  than  that  of  BLX4. 

At  21%  oxygen,  the  ti js  vs.  q"  curves  for  BLX4  and  BLX9  propellants  with 
energetic  plasticizers  present  different  trends  from  those  of  BLX8.  At  the  higher 
energy  fluxes  and  the  higher  oxygen  concentration,  the  delay  times  for  BLX8  are 
significantly  reduced.  However,  for  1  %  and  5  %  oxygen  cases,  the  slopes  are  nearly 
the  same  for  all  BLX-series  propellants.  This  implies  that  at  a  higher  energy  flux, 
the  absorption  of  the  incident  radiative  energy  by  decomposed  gaseous  products 
from  the  GAP  binder  may  become  significant  enough  to  accelerate  reactions  between 
the  product  gases  and  the  oxygen  in  the  ambient  gas. 
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In  contrast  to  those  of  the  BLX-series  propellants,  the  ignition  characteristics  of 
AP/HTPB  composite  propellant  are  nearly  independent  of  the  oxygen  concentration 
in  the  ambient  gas(Fig.l6).  This  is  due  to  the  existence  of  sufficient  oxidizer  in  the 
decomposed  gaseous  products. 

Figure  17  presents  a  more  detailed  description  of  the  ignition  characteristics  of 
BLX8  propellant  at  5%  O2.  For  energy  fluxes  ranging  from  600  to  1000  W j cm2,  two 
different  luminous  flames  appear  at  different  times.  For  example,  at  the  incident 
heat  flux  of  800  W/cm 2  (after  a  period  of  inert  heating),  a  luminous  flame  appears 
briefly,  followed  by  a  dark  period  before  luminous  fiamelets  appear  on  the  sample 
surface.  At  lower  energy  fluxes  (below  600  W/cm2),  no  luminous  flame  distant 
from  the  sample  surface  was  observed.  After  a  relatively  long  exposure  to  the 
laser  beam,  luminous  fl&melets  appear  on  the  sample  surface.  Figure  18  presents  a 
series  of  pictures  recorded  on  the  Spin-Physics  video  camera  at  an  energy  flux  of  700 
W/cm2.  These  photographs  demonstrate  the  existence  of  the  various  zones  depicted 
in  Fig.  17.  The  sample  location,  after  the  appearance  of  luminous  fiamelets  on  the 
sample  surface,  is  indicated  by  two  reticle  lines  in  the  last  picture  of  the  series. 
The  slope  discontinuity  of  the  curve  in  Fig.  17  at  the  critical  energy  flux  (ca. 
600  W/cm2)  is  believed  to  be  caused  by  the  change  of  reaction  mechanism  from 
gas-phase  to  heterogeneous. 

2.4. 1.2  Flame  Observation 

Different  sites  and  shapes  of  flames  were  observed,  depending  on  incident  energy 
flux,  oxygen  concentration  in  the  ambient  gas,  and  composition  of  the  samples. 
Figure  19  shows  the  high-speed  video  pictures  of  luminous  flame  initiation  and 
development  for  a  BLX4  sample  in  a  low  oxygen  concentration.  A  spherical  flame 
which  appears  distant  from  the  sample  surface  expands  to  become  a  peanut-shaped 
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Fig.  17  Ignition  Characteristics  of  BLX8(5%02,latm) 


flame  before  developing  into  a  candle-light  type  flame.  In  contrast  to  this  flame 
development,  in  a  higher  oxygen  environment  a  dome-shaped  luminous  flame  is 
established  close  to  sample  surface.  Figure  20  shows  the  dome-shaped  flame  with 
two  growing  flames  near  the  base  of  the  dome,  which  may  be  due  t,o  the  existence 
of  the  thermocouple  wire.  The  thermocouple  wire  acts  as  a  flame  holder.  After  the 
thermocouple  is  blown  away  from  the  flame,  a  clear  dome-shaped  flame  is  seen  (see 
picture  in  Fig.  20  at  t  =  35  ms). 

During  flame  development,  various  shaped  flames  initiated  from  different  sites 
were  observed,  depending  on  incident  energy  flux,  oxygen  concentration  in  the  am¬ 
bient  gas,  and  composition  of  samples.  The  high-speed  video  pictures  of  ignition  of 
the  BLX4  propellant  sample  (see  Fig.  19)  show  that  at  the  onset  of  light  emission, 
the  luminous  flame  zone  appears  at  a  location  distant  from  the  sample  surface. 
The  flame  expands  into  a  larger  flame,  and  the  lower  boundary  of  the  gaseous  flame 


moves  slowly  toward  the  sample  surface.  Shortly  after  the  appearance  of  the  gaseous 
flame,  the  luminosity  of  the  sample  surface  increases  with  time;  this  indicates  that 
heterogeneous  reactions  occur  on  the  propellant  sample  surface  while  the  products 


of  this  surface  reaction  further  react  in  the  gas  phase.  This  is  drastically  differ¬ 
ent  from  the  flame  development  process  of  BLX8  propellant.  The  most  interesting 
observation  to  be  made  from  these  photographs  is  the  presence  of  localized  bright 
spots  near  the  outer  boundary  of  the  gas-phase  flame  zone  (see  Fig.  19  at  time 
=  43.5,  74.5,  75.0-80.0,  82.5  ms).  These  bright  spots  are  believed  to  be  due  to 
the  entrainment  of  small  RDX  particles  by  the  pyrolyzed  gaseous  species  from  the 
energetic  plasticizer.  Since  BLX4  propellant  contains  a  high  percentage  of  4  fxm 
RDX  particles,  these  particles  can  be  entrained  easily  by  the  pyrolyzed  gases  at  the 
sample  surface.  The  temperatures  of  the  entrained  particles  in  the  dark  zone  are 
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Fig.  18  Flame  Development  of  BLX8  (qn  =  10QW  j cm2 ,b%Oz) 
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Fig.  19  Flame  Development  of  BLX4  ( q "  =  700 W/cm2, 1%C>2) 
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too  low  to  produce  visible  bright  traces.  After  the  particles  have  traveled  a  sufficient 
distance  in  the  gaseous  flame,  they  are  heated  to  an  ignition  condition  and  burn  with 
a  bright  flame.  This  phenomenon  is  not  observed  with  BLX8  propellant,  because 
BLX8  contains  no  energetic  plasticizer  and  because  the  GAP  binder  serves  to  bond 
the  small  RDX  particles  together.  Since  BLX9  has  a  formulation  much  like  that  of 
BLX4  (except  for  plasticizer  type),  the  flame  initiation  and  development  of  the  two 
propellants  are  generally  similar. 

Figure  21  presents  pictures  of  the  flame  development  of  BLX4,  BLX9,  and 
BLX8  propellant  samples  recorded  on  a  conventional  video  camera.  For  all  cases, 
the  oxygen  concentration  in  the  ambient  gas  is  5%.  Since  the  filming  rate  of  the 
conventional  video  is  about  25  pictures  per  second,  detailed  variations  of  flame 
development  could  not  be  captured.  However,  this  recording  does  give  a  broader 
field  of  view  of  flame  development  in  color.  For  the  BLX4  propellant  sample,  the 
luminous  flame  is  initiated  distant  from  the  sample  surface  and  expands  to  reach 
very  near  the  sample  surface.  This  phenomenon  differs  from  the  1%  O2  case  given 
in  Fig.  19.  For  the  BLX9  propellant  sample,  the  flame  is  distant  from  the  sample 
surface  until  the  laser  is  cut  off;  at  this  point  the  flame  disappears.  Fo'  the  BLX8 
propellent  sample,  a  faint  flame  is  established  distant  from  the  sample  surface  and 
the  luminosity  of  the  flame  is  reduced  continously.  Figure  21c  clearly  shows  that  as 
the  gas-phase  flame  weakens,  small  fiamelets  appear  on  the  sample  surface. 

From  the  direct  pictures  recorded  by  the  high-speed  video  camera,  only  the 
luminous  flames  were  observed.  However,  observations  of  the  initial  gasification  and 
movement  of  the  decomposed  gases  are  very  important  in  understanding  the  ignition 
processes.  ’Using  a  schlieren  technique,  the  delay  times  for  the  gas  evolution  and 
the  structures  of  the  primary  flame  are  determined.  Figures  22  show  the  sequences 
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Fig.  22a  Schlieren  Pictures  and  Direct  Pictures 


Wu 


Fig.  22d  Schlieren  Pictures  and  Direct  Pictures 


Fig.  Schlieren  Pictures  and  Direct  Pictures 


Fig.  22h  Schlieren  Pictures  and  Direct  Pictures 
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Fig.  22i  Schlieren  Pictures  and  Direct  Pictures 
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of  gasification  and  the  movement  of  decomposed  gases  recorded  on  the  high-speed 
video  camera. 

2.4.1.3  Thermal  Wave  Structure 

The  surface  temperatures  and  temperature  profiles  in  the  gas  phase  were  mea¬ 
sured  using  the  fine  wire  thermocouples  described  in  the  previous  chapter.  Figure  23 
shows  the  surface  temperature  of  BLX4  in  an  inert  gas  during  radiative  heating  at 
1000  W/cm2.  After  an  initial  transient,  the  surface  temperature  reached  a  constant 
temperature  of  about  480°C.  When  the  laser  is  cut  off,  the  surface  temperature 
decreased  abruptly. 

Temperatures  in  the  gas  phase  were  measured  while  the  surface  of  the  propel¬ 
lant  sample  recessed  during  laser  heating.  The  relative  locations  of  the  thermocou¬ 
ples  to  the  propellant  surface  were  changing  with  time.  Figure  24  shows  the  gas 
temperature  variation  as  the  distance  of  the  thermocouple  from  the  surface  changes 
from  4.7  mm  to  6.5  mm.  As  in  the  surface  temperature  cases,  after  an  initial 
transient  the  gas  temperature  reached  a  nearly  steady  value.  In  the  variation  of  the 
distance  between  the  thermocouple  and  the  receding  surface,  the  spatial  distribution 
of  the  gas  temperature  was  considered  to  be  nearly  uniform.  This  corresponds  with 
the  thermal  wave  structure  in  a  steady-state  combustion.  In  his  study  of  steady- 
state  combustion  of  nitramine  composite  propellants,  Kubota[3]  observed  that  the 
luminous  flame  stood  some  distance  from  the  burning  surface.  He  considered  the 
non-luminous  zone  below  the  luminous  flame  as  a  preparation  zone  similar  to  the 
dark  zone  in  the  flame  of  double-base  propellants.  The  gas  temperature  in  the  dark 
zone  is  nearly  constant. 

Doth  surface  and  gas  temperatures  were  measured  by  means  of  a  thermocouple. 
Before  the  energy  flux  is  applied,  the  thermocouple  wire  is  located  on  the  surface 
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Fig.  25  Surface  and  Gas  Temperature  vs.  Time  (q"  =  lOOOW' / c 
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of  a  propellant  sample.  During  the  initial  phase  of  radiative  heating,  the  thermo¬ 
couple  measures  surface  temperature.  As  heating  continues,  the  surface  of  the  pro¬ 
pellant  sample  recesses  and  the  thermocouple  wire  leaves  the  surface  and  measures 
the  temperature  of  the  gas  phase.  The  location  of  the  thermocouple  wire  during 
the  event  is  monitored  by  the  schlieren  pictures  recorded  on  the  Spin-Physics  high¬ 
speed  video  camera.  Figures  25  and  26  show  the  temperature  traces  for  the  incident 
energy  fluxes  1000  W/cm 2  and  800  W/cm2,  respectively.  After  an  initial  rise,  the 
temperature  reaches  a  constant  value,  about  510°C.  As  heating  continues,  the 
thermocouple  wire  leaves  the  propellant  surface.  Gas  temperature  near  the  surface 
increases  slightly  until  the  second  temperature  rise  reaches  the  preparation-zone 
temperature,  about  1300°C.  It  is  interesting  to  note  that  the  gas  temperature  near 
the  propellant  surface  (about  2  or  3  mm  above  the  surface)  is  nearly  the  same  as  the 
surface  temperature.  This  indicates  that  the  entire  energy  consumed  for  the  surface 
pyrolysis  comes  from  the  incident  laser  energy.  The  contribution  of  conduction  from 
the  gas  phase  is  negligible.  Comparing  the  two  temperature  traces  in  Figs.  24  and 
25,  it  can  be  seen  that  the  surface  and  gas  temperatures  in  the  preparation  zone 
are  nearly  the  same. 

Figure  27  shows  the  temperature  trace  measured  by  a  thermocouple  initially 
located  away  (12  mm)  from  the  propellant  surface.  Once  a  luminous  flame  was 
established  near  the  thermocouple,  the  temperature  increased  abruptly  and  the 
thermocouple  wire  broke.  The  temperature  of  the  luminous  flame  was  too  high  to 
be  measured  by  the  25  n m  Pt-Ptl3%Rh  thermocouple. 


D 


2.4. 1.4  Gas  Analysis 

Gas  samples  from  the  pyrolysis  and  ignition  of  the  nitramine  composite  pro¬ 
pellants  were  analyzed,  using  gas  chromatography  to  determine  the  concentrations 
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Fig.  28  A  Typical  Output  of  Gas  Chromatograph  for  Porapak  Q  and  R  Column 
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of  stable  species.  Various  sizes  of  syringes  were  used  to  take  samples  from  the  test 
chamber.  Since  TVs  is  a  product  gas  from  the  combustion  of  the  nitramine  composite 
propellants,  He  was  used  as  a  carrier  gas  for  the  gas  chromatograph.  Figures  28  and 
29  show  the  typical  outputs  of  the  gas  chromatograph  for  the  test  propellants.  A 
Porapak  Q  in  series  with  a  Porapak  R  column  was  used  to  detect  such  stable  species 
as  Hi,  N2,  CH4,  N20,  NO,  and  CO%.  A  molecular  sieve  (MS  5A)  was  used  to 
detect  O3,  N2,  and  CO.  Many  of  the  peaks  appearing  on  the  outputs  were  identified 
through  calibrations  with  standard  gases  and  cross  examination  with  results  of  a 
GC/MS.  However,  not  all  of  the  peaks  that  appeared  on  the  figures  were  identified. 

The  relative  amounts  of  the  gaseous  species  in  the  product  gases  were  obtained. 
The  weight  percentages  of  the  gaseous  species  were  estimated  by  normalizing  the 
weight  of  the  species  by  the  total  weight  of  the  propellant  sample  consumed  during 
the  event.  As  the  incident  energy  increased,  the  amounts  of  such  stable  species  as 
H2,  C02,  and  CO  increased.  At  a  low  pressure  (ca.  10-3atm),  the  amounts  of 
the  stable  species  decreased  significantly.  Additional  tests  with  the  gas  analysis  are 
necessary  to  understand  the  chemical  reaction  mechanisms  involved  in  the  pyrolysis 
and  ignition  of  nitramine  composite  propellants. 

2.4. 1.5  Sample  Surface  Structure 

Micrographs  of  the  original  and  recovered  BLX4  propellant  samples  are  shown 
in  Fig.  30.  After  a  short  exposure  to  laser  radiation  (before  attaining  light  emssion), 
large  and  small  holes  were  found  on  the  sample  surface  (see  Fig.  30b).  On  the  other 
hand,  surfaces  of  recovered  samples  after  onset  of  light  emission  do  not  exhibit  large 
holes.  Instead,  they  show  more  evenly  distributed  small  holes  surrounded  by  solid¬ 
ified  binder  melt  (see  Fig.  30c).  By  comparing  the  micrograph,  one  can  conclude 
that  during  the  initial  phase  of  laser  heating,  materials  between  RDX  particles  and 


fuel  binder  pyrolyze  and  react  to  create  pressure  buildup  in  the  interface  region. 
Local  pressure  buildup  could  be  the  main  cause  of  RDX  particle  ejection  from  the 
surface.  From  this  point  on,  the  surface  becomes  more  homogeneous. 


2.5  SUMMARY  AND  CONCLUSIONS 


Laser  ignition  of  &  series  of  RDX-based  nitramine  composite  propellants  has 
been  studied  theoretically  and  experimentally.  Observations  and  findings  available 
at  present  are  summarised  as  follows: 

1)  A  comprehensive  theoretical  model  for  radiative  ignition  of  nitramine  composite 
propellants  was  formulated.  A  computer  code  was  developed  to  analyse  the 
theoretical  model. 

2)  A  radiative  ignition  test  setup  using  a  high-power  CO2  laser  was  designed  and 
constructed. 

3)  Ignition  characteristics  of  a  series  of  RDX-based  composite  propellants  were 
studied  experimentally,  using  the  test  setup. 

4)  Delay  time  for  onset  of  light  emission  in  the  gas  phase  is  strongly  dependent 
on  incident  laser  energy  flux,  oxygen  concentration  in  the  ambient  gas,  and  the 
properties  of  propellant  ingredients. 

5)  Propellants  that  contain  energetic  plasticizers  (BLX4  and  BLX9)  are  more 
strongly  dependent  upon  ambient  oxygen  concentrations. 

6)  Critical  energy  fluxes  exist  for  obtaining  distant  gas-phase  flames. 

7)  Particle  ejection  into  the  gas-phase  fame  was  observed;  particle  ejection  is 
believed  to  be  caused  by  local  pressure  buildup  between  RDX  particles  and 
binder  material. 

8)  Micrographs  of  recovered  samples  before  and  after  attaining  light  emission  are 
quite  different.  Large  holes  on  the  propellant  surface  after  tests  indicate  that 
the  RDX  particles  leave  the  sample  surface  before  the  binder  is  pyrolyzed. 

9)  The  effect  of  the  incident  energy  flux  on  the  surface  temperature  of  the  pro¬ 
pellant  samples  is  not  significant  during  pyrolysis  and  ignition. 


10)  A  preparation  Bone  similar  to  the  dark  sone  for  double-base  propellant  flames 
was  observed.  A  very  smalt  temperature  gradient  was  observed  near  the  pro¬ 
pellant  surface  during  ignition.  This  indicates  that  the  energy  consumed  for 
pyrolysis  of  the  condensed  phase  is  supplied  solely  by  the  incident  laser  fluxes. 

11)  As  the  ambient  pressure  and  incident  energy  flux  are  decreased,  the  amount  of 
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III.  CRACK  PROPAGATION  AND  BRANCHING  IN 
BURNING  SOLID  PROPELUNTS 

3. 1  INTRODUCTION 

3.1.1  Motivation  and  Objectives 

Solid  propellant  rookat  safety  is  a  key  faotor  in  several  of 

today's  most  important  oivilian  and  military  propulsion  systems. 

With  both  oivilian  and  military  agencies  presently  pursuing  solid 

rooket  booster  safety  and  ballistic  risk  assesment  programs. 

near-term  advances  in  this  area  are  foreseeable.  Of  particular 

eonoern  is  overpressurization  and  subsequent  catastrophic  failure  of 

a  rooket  motor  due  to  anomalous  burning,  which  may  be  caused  by  the 

presence  of  void3,  cracks,  or  debonded  regions  in  the  propellant 

grain.  Such  regions  are  susceptible  to  further  damage  under 

transient  loading  conditions.  A  theoretical  model  has  been  developed 

to  predict  the  flame  spreading  and  combustion  processes  inside  a 
62 

propellant  cavity  and  various  models  have  also  been  developed  for 
oraok  initiation  and  growth  in  viscoelastic,  inert  mater ial.^”^ 
However,  the  coupling  effects  of  these  two  processes  are  complicated 
and  not  totally  understood.  Based  on  previous  work  by  Kuo  and 
oo-workers®^”^  on  crack  growth  in  burning  propellant  grains,  a  set 
of  governing  dimensionless  parameters  have  been  identified  to  obtain 
empirical  correlations. 

The  primary  objectives  of  this  research  are: 

a)  to  define  a  group  of  meaningful  dimensionless  parameters  for 
characterizing  the  operating  conditions  and  resulting  damage  of  a 
solid  propellant  sample; 
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b)  to  obtain  tost  data  naotssary  for  davaloplng  correlations 
between  the  dlnenalonless  parameters;  and, 
o)  to  establish  correlations,  demonstrating  a  methodology  for 
generating  useful  safety  orlterla. 

3.1.2  Brief  Survey  of  Related  Work 

Early  studies  of  the  oraok  propagation  phenomena  under  burning 

conditions  by  Slefert  and  Kuo  '  provided  good  Insight  into  the  basio 

meohanisms  involved.  They  studied  the  oraok  propagation  process  of  a 

burning  propellant  sample  of  constant  geometry  under  varying 

pressurization  rates.  Experimental  evidence  supported  the  assumption 

that  the  oraok  propagation  velooity  was  limited  by  the  rate  of 

deformation  of  the  oraok  oavlty  walls.  In  this  oase,  a  critioal 

stress  state  was  maintained  in  the  failure  zone  at  the  oraok  tip.  A 

nearly  linear  relationship  between  the  oraok  velocity  and 

pressurization  rate  was  observed, 

V  (m/s)  »  0.086(dp/dt(MPa/s)3°’637  (3-D 

P 

For  low  pressurization  rates  the  oraok  tip  maintained  its  original 

contour.  However,  at  pressurization  rates  on  the  order  of  20  GPa/s, 

the  oraok  tip  appeared  to  square  off  and  a  fan  structure  developed 

suggesting  the  occurrence  of  oraok  branching  and  microstructure 

damage.  It  is  interesting  to  note  that  the  maximum  crack  velocities 

measured  correspond  to  terminal  velocities  of  inert  propellants  and 

70 

rubber  reported  by  Gent  and  Harteny. 

Further  studies  by  Kuo,  Moreci  and  Mantzaras®®*®^  have 
identified  four  different  modes  of  damage  which  were  related  to  the 
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Fig.  31.  Crack  Propagation  Velocity  versus  Pressurization  Rate 
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pressurization  rats  of  tha  oraok  oavlty.  Figure  jl  shows  the 
measured  oraok  propagation  veloolty  as  a  function  of  pressurization 
rate  and  where  eaoh  of  these  modes  was  observed.  Mode  A  represents 
little  or  no  inoreaae  in  oraok  length  and  corresponds  to  very  low 
pressurization  rates.  No  departure  from  normal  burning  is  expeoted 
for  this  oase.  Mode  B  represents  the  propagation  of  a  single  oraok. 
The  rate  of  oraok  growth  Increases  almost  linearly  with 
pressurization  rate  in  the  range  of  1.4  to  15.0  GPa/s.  Although 
burning  is  enhanoed  at  the  oraok  tip,  the  overall  inorease  in  burning 
due  to  the  damage  is  not  significant.  Mode  C  represents  an  Initial 
single  oraok  propagation  with  looal  branching  at  various  locations. 
This  was  observed  for  Intermediate  pressurization  rates  and  has  the 
potential  of  significantly  increasing  the  burning  surface  wreu.  Mode 
D  represents  oraok  branohlng  from  the  Initial  crack  tip  looation  in 
either  descrete  raacrocracks  or  a  fan  of  numerous  smaller  cracks. 

This  type  of  damage  was  characteristic  of  higher  pressurization  rates 
and  also  can  significantly  increase  the  burning  surface  area.  Some 
of  the  data  from  these  experiments,  as  well  as  newly  acquired  data, 
were  used  in  the  ourrent  analysis. 
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3.2  EXPERIMENTAL  SETUP 


3.2.1  Seleotlon  of  Sample  Propellant  Type 

In  selecting  a  propellant  sample  for  this  research,  the  following 
orlterta  were  considered: 

a)  the  propellant  must  be  available  for  a  basic  study  in  a 
university  environment; 

b)  Its  oombustion  behavior  must  be  well  oharaoterlzed; 

o)  its  mechanical  properties  and  failur«  criterion  under  non-bur ing 
conditions  must  be  known;  and, 

d)  propellants  of  the  same  family  must  be  available  so  that  similar 
tests  can  be  conducted  to  determine  the  effect  of  burnig  rate  and 
oxidizer  site. 

Two  composite  propellants  (ARC  5051  and  ARC  4525)  were  selected.  Each 
of  these  contain  271  HTPB  and  731  AP  with  average  oxidizer  partlole 
sizes  of  200Um  (ARC  5051)  and  20Um  (ARC  4525).  Due  to  the  similarities 
in  composition  thermodynamic  propeties  such  as  enthalpy,  density, 
burning  rate,  etc.,  are  nearly  equal.  Their  combustion  behavior  and 
mechanioal  properties  have  been  well  characterized.  In  addition, 
propellants  of  similar  compostion  but  with  different  burning  rates  are 
available  for  future  studies.  Although  these  propellants  have  a  lower 
solids  loading  than  the  typical  high-energy  solid  propellants  used  in 
most  of  todays  advanced  rocket  propulsion  systems,  once  the  basic  method 
for  determining  oraok  oombustion  correlations  has  been  established, 
similar  studies  oan  be  made  for  other  propellants. 


3.2.2  Sample  Configuration 


The  two-dimensional  sample  configuration  used  in  previous  studies 
is  adopted  for  this  study  (?  ;g ,  32).  This  general  geometry  is  similar 
to  one  p<'int  of  a  star  grain  or  to  a  segment  of  a  grain  containing  a 
void.  This  configuration  provides  an  initial  crack  of  variable  length 
and  propagation  along  the  horizontal  axis.  The  geometric  parameters  iJj 
and  i|>  can  also  be  varied  to  represent  different  local  boundary 
conitions  cf  the  damaged  zone. 

3.2.3  Test  Setup 

The  test  setup  was  designed  to  induce  combustion  and  mechanical 
damage  in  a  localized  zone  of  a  pre-cut  propellant  sample  by  subjecting 
it  to  high-pressure,  high-temperature  gases  and  to  interupt  the  buring 
by  rapid  depressurization  for  sample  recovery.  In  these  tests 
variations  of  several  independent  parameters  were  made,  including 
pressurization  rate,  initial  crack  length,  sample  geometry  and 
propellant  composition.  The  damage  was  then  measured  in  terms  of  cracw 
propagation  velocity,  extended  crack  length  and  number  of  cracks. 

3.2.3. 1  Test  Chamber 

The  test  rig  (Fig  33)  consisted  of  a  high-pressure,  windowed 
chamber  with  the  following  subsystems; 

a)  a  driving  motor  to  pressurize  the  free  volume  of  the  chamber; 

b)  a  rupture  disk  and  venting  port  to  rapidly  depressurize  the 
chamber;  and, 

c)  a  ni<-"ogen  injection  system  to  cool  the  exposed  surfaces  and 
reduce  the  cnance  of  sample  reignition. 
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Rupture  Disk 


Fig.  33.  Test  Chamber 
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3. 2.3.1. 1  Driving  Motor 

The  driving  motor  contained  four  key  components  (Fig  34).  The 
ignition  primer  (type  M52A3B1),  when  supplied  with  three  to  four  volts 
DC,  would  fire  and  in  turn  ignite  the  propellant  shavings.  The  shavings 
burn  extremely  fast  due  to  their  very  high  burning  surface  area.  Their 
purpose  was  to  supply  hot  gases  and  boost  the  pressure  in  the  segment 
containing  the  main  igniter  charge.  The  igniter  charge  would  then  burn, 
exhausting  through  the  perforated  nozzle  into  the  main  chamber.  By 
varying  the  mass  of  the  main  igniter  charge,  pressurization  rates 
ranging  from  near  zero  to  50  Gpa/s  were  achieved. 

3.2.3. 1.2  Rupture  Disk  and  Venting  Port 

When  the  chamber  pressure  exceeds  a  predetermined  value  the  brass 

rupture  disk  bursts,  allowing  the  high  pressure  gases  to  escape  through 
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the  venting  port.  As  described  by  De  Luca,  dynamic  extinction  can  be 

achieved  by  exceeding  a  critical  depressurization  rate.  This  rate 

increases  as  chamber  pressure  increases  and  is  also  substantially 

affected  by  propellant  composition.  Due  to  its  complexity,  no  attempt 

was  made  to  characterize  the  extinction  properties  of  the  test 

2 

propellants.  Instead,  an  exit  port  area  of  2.54  cm  ,  based  on 
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experimental  data  of  Kuo  and  Moreci,  used  with  rupture  pressures 
between  7.0  and  15.0  MPa  proved  to  be  sufficient. 

3.2.3. 1.3  Nitrogen  In.iection  System 

After  the  rupture  disk  bursts  and  the  chamber  pressure  is  reduced, 
hot  spots  on  the  chamber  walls  could  reignite  the  propellant  surfaces 
and  the  sample  would  be  completely  consumed.  In  an  effort  to  improve 
the  reliability  of  interrupting  the  burning  process,  nitrogen  was 


Injected  following  depressurization  to  oool  the  surfaces  of  the  free 
volume. 

3. 2, 3. 2  Data  Acquisition  System 

The  data  needed  from  these  tests  to  calculate  key  parameters 
consisted  of  a  pressure-time  trace  of  the  craok  cavity  as  well  as  a 
visual  record  of  the  crack  itself  as  it  propagated.  A  complete 
schematic  of  the  data  acquisition  system  is  shown  in  Fig  35. 

3. 2. 3. 2.1  Pressure  Measurements 

The  pressure  measuring  device  selected  was  a  miniature 
piezoelectric  quartz  transducer  (Kistler  Model  601B1).  This  transducer 
has  a  10%  to  90%  rise  time  of  3.0  microseconds  over  the  range  of  0  to 
15.000  psi.  This  was  more  than  adequate  for  measuring  pressurization 
rates  of  up  to  50  GPa/s  with  rise  times  on  the  order  of  1  millisecond. 

To  protect  the  transducers  from  the  high  temperature  product  gases,  they 
were  recess-mounted  in  water-cooled  adapters.  An  insulated 
high-impedence  cable  with  a  nominal  capacitance  of  30  pF/ft  was  used  to 
transmit  the  signal  from  the  transducer  to  the  charge  amplifiers.  The 
signal  attenuation  under  the  test  concH  *■  ns  was  less  than  one  percent. 

The  charge  amplifiers  (Kistler  Model  SOME)  output  the  transducer 
charge  signal  au  an  amplified  voltage.  Calibration  is  achieved  in  the 
charge  amplifier  by  supplying  a  DC  voltage  of  equal  magnitude  to  that  of 
the  pressure  transducer  sensitivity.  The  output  voltage  of  the  charge 
amplifier  is  then  adjusted  to  satisfy  the  equation 


-V.  x  lOOOpF 

v _ IE _  . 

out  "  Range  x  Sensitivity 


(3-2) 


A  typical  range  setting  of  2000  psi/volt  resulted  in  a  calibration 
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output  of  0.5  volts.  A  long  time  constant  was  used  to  provide  a 
continuous  signal. 

The  output  from  the  charge  amplifiers  was  recorded  on  a  Physical 
Data  Model  515A  Transient  Waveform  Recorder.  It  reoorded  the  voltage 
signals  at  sampling  intervals  of  5  to  50  microseoonds  for  the  initial 
pressurization  period  and  then  shifted  to  a  sampling  rate  of  2.5  to  5 
milliseconds  in  order  to  extend  the  total  time  reoorded.  The  limiting 
factor  was  its  maximum  oapaoity  of  4096  data  points.  The  recording 
mechanism  was  triggered  by  a  small  signal  from  the  transducer  in  the 
craok  cavity.  A  negative  delay  of  512  data  points  was  used  to  insure 
the  beginning  of  the  event  was  recorded.  The  recorded  signal  could  then 
be  played  back  to  a  Nicolet  Digital  Osoillosoop:*  for  analysis  and 
recording  on  floppy  disk  or  to  an  X-Y  plotter  for  a  hardcopy. 

3. 2. 3. 2. 2  High  Speed  Film  Recording 

A  Hycara,  16  mm  camera  was  used  to  record  the  ignition,  combustion, 
and  crack  propagation  of  the  oropellant  sample  through  a  veiwing  port  in 
the  chamber.  Typical  framing  rates  were  5,000  frames  per  second  in  a 
half  frame  format  resulting  in  10,000  pictures  per  second.  In  order  to 
insure  the  event  occured  before  the  roll  of  film  was  completely  exposed, 
the  camera  was  used  to  trigger  the  ignition  of  the  driving  motor.  The 
camera  was  started  and  when  a  prescribed  footage  of  film  was  exposed, 
allowing  the  camera  to  reach  the  set  film  speed,  an  internal  switch 
would  close  the  igniter  circuit.  This  would  also  trigger  a  common  time 
pulse  which  was  recorded  both  on  the  film  and  on  the  Physical  Data 
Recorder.  Another  light  pulse  with  a  frequency  of  1  KHz  was  also 
recorded  on  the  film,  to  one  side  of  the  main  image,  in  order  to  verify 
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the  framing  rat* 


3.2.4  Test  Procedure 

Appendix  F  contains  a  complete  oheckllst  and  sample  data  sheet  used 
for  the  lnterupted  craok  propagation  tests.  The  following  is  a  brief 
description  of  the  key  steps  in  a  typioal  test  firing: 

a)  measure  the  main  Ignitor  charge  and  shavings  and  assemble 
ignitor.  Measure  the  resistance  of  the  primer  and  short  the  primer 
oircult; 

b)  measure  and  record  the  craok  geometry  of  the  propellant  sample. 
Cover  all  surfaces,  except  crack  cavity,  with  flame  retardant  grease  to 
prevent  flame  spreading  over  those  surfaces.  Install  the  sample  and 
bolt  on  the  window  and  retainer; 

c)  Install  pressure  transducers  and  calibrate  the  charge 
amplifiers.  Connect  amplifier  outputs  to. the  Physical  Data  recorder; 

d)  set  up  the  high-speed  camera  and  test  ingition  circuit; 

e)  put  charge  amplifiers  in  the  OPERATE  position,  reset  the 
Physical  Data  recorder  and  start  the  camera; 

f)  when  the  rupture  disk  bursts,  open  the  nitrogen  line; 

g)  turn  off  charge  amplifiers,  store  and  plot  the  pressure  traces 
and  disassemble  the  chamber  to  examine  the  recovered  sample. 


3.3  THEORETICAL  ANALYSIS 


3.3.1  Crack  Sample  Geometry 

Figure  36  shows  the  damaged  sample  geometry  considered  in  this 

analysis.  From  this  geometry,  some  key  parameters  oen  be  defined. 

First,  and  4^  are  Irect  geometric  factors  af feotlng  the  oraok 

propagation  prooess,  and  the  number  of  oraoks  n,  is  a  key  factor  in 
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defining  the  degree  of  damage.  The  affected  volume,  or  damaged-zone 
volume  is  defined  as 


Ya  «  irL^b  (3-3) 

The  burning  surfaoe  surface  area  in  the  affeoted  region  is  given  by 
Ajj  a  2Tr(n+1)LDb  (3-4) 

For  this  analysis,  a  specific  burring  surface  area  is  defined  as 

A  _  *b  .  2(n+1)  (3-5) 


The  functional  forms  of  equations  (3-3)  and  (3-5)  are 

Va  *  VLD*  b)  <3-6 ) 

As  a  As(Ld>  n)  (3-7) 

3.3.2  Determination  of  Mechanical  Properties 

A  series  of  mechanical  tests  was  conducted  by  Harbert  and 
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Schapc'y  to  determine  the  mechanical  properties  of  the  two  propellants 
(ARC  5051  and  ARC  4525)  under  non-burning  conditions. 
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3. 3. 2.1  Relaxation  Modulus 


-Jniaxlal  constant  strain  tasts  on  propellant  bar  specimens  were 
par fora ad  at  various  taaparaturas.  Relaxation  moduli  varaua  tlma  data 
ara  plottad  for  aach  hamper atura.  In  order  to  obtain  a  continuous  ourva 
for  all  taaparaturas,  the  relaxation  moduli  ara  plottad  versus  a  shifted 
time  £,  where 

K  «  t/ay  (3-8) 

Tha  tlma  shift  faetor  a^  was  datarainad  based  on  a  reference  temperature 
of  297  K  (l.e.,  aT  •  1.0  at  297  K).  Figure  37  shows  tha  funotlonal 
relationship  between  temperature  and  tlma  shift  faotor.  The  resulting 
master  relaxation  curve  for  eaoh  propellant  is  shown  in  Figure  38. 

3.3. 2.2  Crack  Propagation  Velooltv 

Craok  growth  tests  were  conducted  on  strip  specimens  with  a  preout 
crack  (Fig.  39)  at  different  orosshead  rates.  The  oraok  velooity  and 
strain  were  measured  directly  and  the  stress  intensity  was  calculated 
using 

Kj  «  o  h/2(1-v2)  (3-9) 

where  Poisson's  ratio  v  *  0,5  and  o  is  the  stress  far  from  the  crack  tip 
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exerted  by  the  crosshead.  Figure  40  shows  the  relationship  between 

crack  propagation  velocity  and  stress  intensity  faotor  for  ARC  4525 

propellant.  It  oan  be  seen  that  the  lines  of  constant  strain  are 

parallel  (l.e  independent  of  strain).  This  faot  suggests  the 

applicability  of  linear  viscoelastic  theory  to  characterize  the  craok 
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propagation  velooity.  The  results  of  this  theory  can  be  expressed  by 
the  following  power  law 
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Fig.  39.  Test  Configuration  for  Crack  Velocity  Characterization 


(3-10) 


In  this  equation  q  is  the  strain-independent  slope  obtained  from  Figure 
40,  and  C  is  a  material  parameter  which  varies  with  strain.  The 
functional  representation  of  equation  (3-10)  is  then 

vp  =  VKi’  e.  «T)  <3-10 

3. 3.2.3  Fracture  Energy  and  Related  Parameters 

From  an  idealized  crack  propagation  model  for  viscoelastic  media  as 
used  by  Knauss^3  and  Sehapery^”^,  the  fracture  energy  (r)  can  be 
defined  as  the  energy  required  to  elongate  a  segment  of  area  dA  to  the 
point  of  failure.  This  quantity  has  the  functional  relationship 

T  *  r«r  Vp,  q>fin)  (3-12) 

where  cr^  is  the  maximum  stress  in  the  failure  zone.  The  stress 
intensity  factor  is  generally  a  function  of  the  geometry  of  the 
continuum  in  which  the  failure  zone  is  embedded  and  the  stress  in  that 
zone,  represented  by  chamber  pressure  in  this  application.  Thus 

KI  s  KI(ao’  V  ^i*  P)  (3-13) 

and 

K.  =  Kt  (3-14) 

iC  A 

76 

for  crack  initiation.  From  Langlois  and  Gonard,  the  extended  crack 
length  can  be  related  to  stress  intensity  factor  and  the  maximum  stress 
in  the  failure  zone 

LD  ■  lD«I>  <3-,5> 
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3.3.3  Energy  Balance  on  the  Control  Volume 


An  energy  balanoe  is  applied  to  the  control  volume  (Figure  41 )  in 
order  to  relate  key  parameters  in  the  oombustion  process.  The  basic 
assumptions  are: 

a)  The  heat  loss  from  the  control  volume  to  the  surroundings  is 
negligible  since  the  damage  occurs  in  an  extremely  short  time  interval 
(of  the  order  of  10  ms). 

b)  The  combustion  gases  of  the  driving  motor  are  considered  to 
follow  the  porfect  gas  law.  This  assumption  is  valid  since  the  gas 
temperature  is  high  (2000  K)  and  the  maximum  chamber  pressure  is  less 
than  35  MPa. 

c)  mQut  =  0  for  the  time  of  interest  since  the  initial  damage 
occurs  before  the  rupture  disk  bursts. 

d)  The  gas  phase  control  volume  can  be  treated  as  a  constant  since 
the  volume  increase  due  to  flow  work  and  surface  burning  is  extremely 
small  compared  to  the  total  control  volume  during  the  time  of  interest. 
The  conservation  of  energy  for  the  gas  phase  can  be  written  as 


*  &  h_  -  A  .  c  T  -  0  -  P 


3t  “  ign  t,ign  “p‘  p  out^p 


av 

_ c 

a  t 


(3-16) 


Using  the  basic  assumptions,  equation  (3-16)  reduces  to 


yrl  3t  "  ^ign^t.ign  +  ^p^p 


(3-17) 


After  replacing  according  to 


np  "  prb^bt 


(3-18) 


where  the  total  burning  surfaoe  area  oan  further  be  related  to  the 
orack  geometry 

St 3  St(V  S'  S'  S'  V  (3-19) 

Equation  (3-17)  then  takes  on  the  functional  form 


3P  3P 


it  ■  it‘V»-  V'  S.ign*  S'  S'  S'  S'  Pp.  S'  S'  V  (3-20) 

The  burning  rate  under  a  pressure  transient  has  been  characterized 


by  Krier77  as 


na 


r.  (t)  =  r  [1  +  * - B-r 


(3-21) 


P(t)r  dt 
o 


where  rQ  is  the  steady  state  burning  rate  according  to  the  Saint-Robert 
78 

law,  oip  is  the  solid  phase  thermal  diffusivity,  and  n  and  $  are 

numerical  constants.  Since  the  thermal  diffusivities  were  the  same  for 
73 

both  propellants  it  also  was  treated  as  a  numerical  constant.  The 
instantaneous  burning  rate  can  therefore  be  related  as 

rb  a  rb(P*  3P/3t)  (3-22) 

Since  the  flow  from  the  driving  motor  is  choked  and  the  ignitor 
charge  is  made  of  the  same  propellant  as  the  sample,  the  propellant 
enthalpy  can  be  related  to  the  stagnation  enthalpy  of  the  igniter  gases 
as 


h 


P 


2  . 

Y~1nt,ign 


(3-23) 


The  propellants  are  essentially  incompressible  and  of  equal 


density,  therefore 


constant 


(3-24) 
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From  evidence  obtained  in  oraok  propagation  experiments,  crack 
propagation  velocity  in  the  burning  solid  propellants  is  related  to  the 
following  parameters: 


VP  *  vp^p/^*  p*  mign’  r*  V  V  Aa,  aQ,  aT) 


(3-25) 


3.3.4  Dimensional  Analysis 

Equations  (3-6),  (3-7),  (3-1 1  >-(3-15),  (3-20)  and  (3-22M3-25) 
relate  the  twenty-three  parameters  (y.was  assumed  to  be  a  numerical 
constant  for  the  combustion  product  gases)  affecting  the  crack 
propagation  process.  These  twelve  equations  may  be  reduced  to  one 
equation  involving  twelve  of  the  more  significant  parameters,  namely 

vp  *  Vp<»P/3t.  V  V  ly  Kjc.  Va>  *s,  Ljj.  *t>  a,)  (3-26) 

3.3.4. 1  Arrangement  into  Meaningful  Dimensionless  Groups 
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The  Buckingham  Pi  Theorem  indicates  that  there  exist  nine 
independent  dimensionless  parameters  involved  in  the  crack  propagation 
process.  Three  of  the  key  parameters  are  already  dimensionless,  thus 
i|>o,  and  were  chosen  as  dimensionless  parameters  describing  the 
geometry  of  the  continuum  around  the  damaged  zone  and  the  effect  of 
initial  propellant  temperature  on  mechanical  properties.  In  order  to 
develop  additional  meaningful  groups,  the  parameters  which  desoribe  the 
degree  of  damage  (i.e.,  L^,  Va,  Ag  and  V  )  were  examined  individually. 

A  desireable  dimensionless  parameter  defining  the  extended  crack 


length  was  the  percentage  increase  in  overall  crack  length.  This 
dimensionless  parameter  was  defined  as 


(3-27) 


it  _  _  Length  of  extended  maorooraok 

A  “  aQ  Length  of  the  Initial  eraok 

From  the  parameters  V  and  A  ,  two  meaningful  dimensionless  groups 

8  S 

oould  be  formed, 

V  A  ^ 

^  _  a  s  _  Rate  of  increase  of  oraok  surface  area  (3-28) 

B  "  ~  Rat®  of  increases  of  cross-sectional 

UD  *  area  of  the  affected  region 

ftg  «  LDAS  -Ian,  number  of  macrocracks  of  length  LQ  (3-29) 
The  crack  propagation  velocity  was  normalized  to  the  normal 
regression  rate  of  the  surface  due  to  burning, 

V 

u  _ Crack  propagation  Velocity _  (3-30) 

6  “  r0  Burning  rate  of  the  main  propellant 


The  remaining  dimensionless  parameters  relate  the  input  energy 
from  the  propellant  and  the  presurization  rate  to  the  critical  values 
for  crack  propagation  and  were  grouped  as  follows: 

h  o  r  A  irL^  Rate  dermal  ®n®rgy  input 

^  _  pPprbAs7rLd _ from  propellant  product  gases 

E  ”  V  [K^  (1-v^)/E  ]  ~  s"ra*n  ®nergy  release  rate 

pL  Ic  ~  refJ  required  for  crack  propagation 

dP  ^  3/2  Rate  of  change  of  pressure  force  exerted 

u  _  dt  a  _  on  the  crack  surface _ 

Kj  rfa  “  (critical  stress  for  fracture) x( rate  of 

C  increase  of  flow  area  in  the  damaged  zone) 

The  additional  constants  such  as  density,  sample  thickness  and 
Poisson's  ratio  were  added  to  give  the  groups  their  respective  physical 
mearing. 

3. 3. 4. 2  Functional  Relationships  Between  Dimensionless  Parameters 
The  Dimensionless  parameters  irA,  tt0,  irc,  irD  characterize  the 
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(3-31 ) 


(3-32 ) 


degree  of  damage  in  the  propellant  sample.  The  driving  force  oausing 
this  damage  (l.e.  the  Initial  pressurization  rate  and  thermal  energy 
input  from  the  propellent  gases)  are  oharaoterized  by  the  dimensionless 
parameters  and  These  parameters  were  used  together  with  ^0,  ^ 
and  aT  for  constructing  the  following  emperioal  correlations  to  predict 
the  oraok  propagation  and  branohing  processes: 

UA  *  TTp,  1^0 » 

"b  *  1tb<1tE '  V  *»•  +f  *T>  (3.33) 

*  ^C^E*  ^F*  *o*  ^i*  aT^ 

ttd  =  V  V  V  aTJ 
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3.4  RESULTS 


3.4.1  Determination  of  Dimensionless  Parameters 

Several  tests  were  oonduoted  using  the  eparatus  previously 
desorlbed.  From  the  pressure  data,  film  data  and  known  meohanioal  and 
thermodynamio  properties,  values  for  eaoh  of  the  parameters  oould  be 
deduced. 

3. 4. 1 . 1  Pressure  Loading  in  the  Damaged  Zone 

Figure  42  shows  the  pressure-time  trace  at  three  key  looations  in 
a  typioal  test.  Looation  1  corresponds  to  the  chamber  free-volume  and 
looatlon  2  corresponds  to  the  initial  orack  tip  location  and  the  3rd  is 
further  along  the  craok  propagation  direction.  The  first  curve  is  the 
driving  pressurization  rate  used  in  calculating  the  dimensionless 
parameters  since  it  is  generally  independent  of  the  Initial  orack 
configuration.  The  pressure  trace  for  the  second  location  illustrates 
the  increased  pressurization  rate  end  peak  pressure  experienced  at  the 
tip  of  an  opening  crack.  This  may  be  caused  by  the  delayed  opening  of 
the  crack  and/or  the  increased  burning  rate  in  the  failure  zone.  As 
seen,  the  peak  pressure  in  this  case  was  13.3  HPa  (2000  psi)  when  the 
rupture  disk  burst.  The  chamber  pressure  then  dropped  rapidly,  and  the 
sample  was  recovered.  Figure  43  is  a  close-up  photograph  of  the 
damaged  zone  showing  the  bracnching  which  occurred. 

Figure  44  is  a  similar  pressure  trace  for  another  test.  The 
rupture  disk  burst  earlier  in  this  case  (7.0  MPa),  whioh  r.ay  have 
contributed  to  the  reignition  and  burn-out  of  the  sample.  The  early 


Fig.  42.  Pressure-Time  Trace  for  DNICP  27 
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rupture  appeared  to  affeot  the  peak  prtaaura  at  the  crack  tip  aa  it 
ooQurrtd  during  the  depreaaurlxatlon  of  the  free  volume. 

3.0.2  Interpretation  of  the  Film  Data 

The  degree  of  damage  to  the  propellant  aample  was  eharaoterized  by 
exaailnlng  the  film  reoord  of  the  event,  A  typloal  film  record  of  a 
Node  B,  single  oraok  propagation  la  shown  in  Figure  *5.  Table  IV 
oontalns  a  frame-by-frame  deoription  of  the  event.  The  instantaneous 
bright  flame  oontoura  caused  by  burning  in  the  crack  are  used  to  define 
the  oraok  shape  and  length s  From  the  oontour  (Figure  46),  the  oraok 
tip  looation  was  measured  as  a  function  of  time.  Thus  the  oraok 
velocity,  extended  length  and  number  of  maorooraoks  were  deduced.  In 
this  test  the  initial  oraok  length  was  longer  than  in  previous  tosts. 
Figure  47  shows  that  the  luminous  flame  front  moved  through  the 
initially  olosed  region  and  the  oraok  propagation  region  at  the  same 
rate.  The  sharp  reduction  in  propagation  rate  was  caused  by  the 
depressurlxation  of  the  ohamber.  A  typloal  instantaneous  bright  flame 
oontour  for  a  branching  experiment  is  shown  in  Figure  48. 

3.4. 1.3  Determining  Other  Key  Parameters 

In  order  to  calculate  each  dimensionless  parameter  the  values  of 

lb»  KI  *nd  Eref  rauat  be  deternin®<1*  To  represent  the  normal 
C 

regression  of  the  propellant  surface  in  relation  to  the  oraok 
propagation,  a  steady  state  burning  rate  was  calculated  using  Saint 
Robert’s  law^® 

rb  =  aPn  (3-34) 

where  P  was  the  pressure  at  the  time  of  oraok  initiation  and  a  and  n 
are  burning  rate  constants  given  in  Table  V.®J 


Table  IV.  Film  Interpretation  of  Crack  Propagation  Process 
Picture  Ho.  Time  (ms) _ Event _ 

1-6  0-0.6  Igniter  gases  fill  crack  cavity. 

7-12  0.6-1. 2  Flame  penetrates  precut  crack. 

13-16  1.2-1. 6  A  single  crack  propagates.  High 

pressure  gas  continues  to  flow 
into  crack  region. 

17-28  1.6-2. 8  Crack  cavity  is  clearly  seen. 

Propagation  has  stopped. 

29-31*  2. 8-3. 4  Rupture  disk  bursts  and  flame 

intensity  is  reduced. 


35 


3.4- 


Propellant  sample  is  not  extinguished' 
and  burns  completely. 


■V.V.'v 


3P/3t  -  12.0  GPa/s 

V  -  35. A  b/s 
P 

n  ■  1 


Fig.  46.  Instantaneous  Bright  Film  Contour  From  DNICP  40 
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Table  V.  Steady  State  Burning  Rate  Constants 


Using  the  data  obtained  from  numerous  experiments,  several  sets  of 
values  for  the  nine  dimensionless  parameters  have  been  deduced  (Table 
VI),  With  this  data,  the  functional  relationships  of  these  parameters 
were  obtained  statistically  using  a  multiple  regression  analysis 
program.  The  results  can  be  expressed  as 


IT,  .  20.697  ,  0.107  0.181  I0'6"5 

A  c.  r  l  1  o 


(3-36) 


uB  =  6.16X10"7  7TE°-0657rF0.254aT-0.78(sin^)-6.il5(sin^r0.985 


(3-37) 


»c  =  7.389  uE0-27V'04V,021(si^i)1,11'(aln^o)"0*0n5 


(3-38) 


*D  =  10.9x10^  wE0-2;v*07S0*1,78(3i^i)2*05(si^or0'117 


(3-39) 


It  ic  important  to  note  that  although  the  exponents  of  and 
are  small  compared  to  the  other  exponents,  the  values  of  these  two 


parameters  are,  in  general,  large  compared  to  sin^,  sin^Q  and  a^,  and 

therefore  have  a  dominant  effect. 

In  order  to  oheok  the  reliability  of  these  predictions,  equations 

(3-36)  through  (3-39)  are  plotted  with  the  values  obtained  from  test 

data  (Figures  49  through  52).  These  figures  show  that  the  predictions 

adequately  oharaoterize  the  experimental  data.  The  multiple  regression 

2 

analysis  calculated  R  values  greater  than  89%  for  all  correlations. 

2 

R  is  defined  as  the  sample  multiple  coefficient  of  determination  given 
by 


R 


2 


S«1-  v2 

z(V  r)2 


(3-40) 


where  are  the  measured  values,  T  is  the  mean  value  of  and  Y^  is 
the  predicted  value  of  Y^ 


3.4.3  Disoussion  of  Results 

The  expressions  for  •nrA,  fig,  ir  and  characterize  the  damage  in 

the  sample  in  terms  of  the  extended  macrocrack  length,  increase  in 

burning  su  .ce  area,  number  of  macrocracks,  and  crack  propagation 

velocity.  From  equations  (3-36)  through  (3-39),  the  effects  of  the 

operating  conditions  characterized  by  ljj  and  a^  on  the 

sample  damage  can  be  analyzed  individually.  The  equations  show  that 

all  damage  parameters  increase  with  ffg  and  Up.  The  pressurization  rate 

OTp)  is  par^culai*  lominant  in  extending  the  crack  length  and 

increasing  the  specific  surface  area.  This  agrees  with  previous 
67-69 

studies  which  showed  an  increase  in  damage  for  higher  rates  of 
thermal  energy  input  higher  pressurization  rates. 

Changes  in  cnange  the  confinement  of  the  damaged  zone  allowing 
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Fig.  4S.  Data  Correlation  for  the  Extended  Crack  Length  IIx 
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Fig.  52.  Data  Correlation  for  the  Crack  Propagation  Velocity  JId 
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varying  degrees  of  sample  dafornatlon.  The  aquations  predlot  that  a 
graatar  anglt  will  inoraaaa  oraok  growth  while  aaoh  of  tha  othar 
parameters  will  daoraaaa.  Undar  thasa  conditions  tha  angla  at  tha 
oraok  tip  is  likaly  to  ba  graatar  raaulting  in  a  highar  looal  stress 
concentration.  Thus,  tha  single  oraok  propagation  node  is  more 
probable.  With  a  longer,  single  oraok,  equation  (3-37)  indioates  that 
thasa  two  affaots  oontrlbute  to  a  decrease  in  ffg.  Also,  higher  sample 
dafornatlon  work  implies  lass  energy  is  available  to  propagate  tha 
oraok,  reduoing  Vp. 

The  affaots  of  initial  propellant  taaperature  on  the  oraok 
propagation  prooess  can  be  vary  important.  At  lower  temperatures 
(higher  values  of  a^)  tha  propellant  becomes  more  brittle,  lowering  its 
fracture  energy.  Higher  oraok  velocities  and  longer  oraoks  are 
expeoted.  However,  the  range  of  initial  temperatures  in  this  study  was 
limited  by  the  test  conditions.  Broader  ranges  are  needed  to  more 
aoourately  predict  its  effect.  ^  was  sssentially  unchanged  In  these 
tests  but  was  inoluded  for  the  completeness  of  the  correlations. 

The  agreement  of  these  predictions  with  previous  observations  and 
basio  principles  in  fracture  mechanics  indicates  that  the  methodology 
for  their  development  is  valid.  Similar  correlations  for  different 
propellant  types  with  different  oraok  configurations  oould  be  developed 
using  this  method  to  determine  guidelines  for  the  safe  operation  of 
high  performance  rooket  motors. 
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3.5  SUMMARY  AND  CONCLUSIONS 


Two  types  of  con post to  propollonts  wort  solootod  for  this  study 
(ARC  5051  and  ARC  4525).  Thoy  havt  similar  composition  but  dlfforont 
oxidizor  partiolo  sizts  and  moohanioal  properties.  Several  tests  wore 
oonduotod  under  different  loading  conditions  and  sapmle  geometries. 

The  damage  to  the  sample  was  measred  in  terms  of  oraok  length,  number 
of  oracks  and  oraok  propagation  velooity.  Using  a  dimensional  analysis 
of  the  key  variables  Involved  in  the  oraok  propagation  phenomena,  a  set 
of  governing  dimensionless  parameters  were  established.  These 
parameters  oan  provide  some  Insight  into  the  cause  and  effect 
relationship  between  the  momentum  and  energy  input  to  a  propellant 
sample  with  existing  flaws  and  the  resultant  damage. 

Correlations  of  the  dimensionless  parameters  were  obtained  by 
performing  a  multiple  regression  analysis  on  the  test  data.  The 
driving  parameters,  and  Up  have  positive  relationships  with  the 
damage  parameters  Ug,  ard  u^.  This  indicates  that  higher 
momentum  and  energy  inputs  to  the  oraok  cavity  will  cause  a  greater 
degree  of  damage.  The  geometry  and  initial  temperature  of  the 
propellant  sample  have  additional  effects  characterized  by  the 
correlations.  These  corelations  are  only  valid  for  the  small  ranges  of 
temperature  and  crack  cavity  angles  studied,  however,  additional  tests 
could  expand  the  applicability  and  validity  of  this  method.  Enough 
data  has  been  obtained  to  predict  the  strong  dependence  on  the  rate  of 
energy  input  and  pressurization  on  orack  propagation.  The  good 
agreement  with  test  data  and  reasonable  physical  interpretation 
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suggests  that  the  methodology  outlined  in  this  study  can  provide  a 
basis  for  future  development  of  similar  correlations  defining  safe 
operating  conditions  for  solid  propellant  rockets. 

Future  research  in  this  area  should  include: 

a)  development  of  correlations  for  high-energy  propellants  used  in 
actual  propulsion  systems  today,  including  high-elongaticn  propellants; 

b)  extension  of  the  present  analysis  to  include  propagation  and 
flame  spreading  along  a  motor  casing  due  to  an  initially  debonded  area, 
and; 

c)  integration  of  a  crack/ debond  model  into  the  Rocket  Performance 
Prediction  Code. 
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Appendix  A 


MEASUREMENT  OF  CO.  LASER  ENERGY  FLUXES 

It  is  very  important  to  define  the  energy  fluxes  delivered  by  the  incident  laser 
beam  for  ignition  studio*.  In  order  to  define  the  COj  laser  energy  fluxes,  it  is 
necessary  to  measure  the  laser  intensity  distribution  precisely.  In  this  study,  three 
different  devices  were  used  for  the  measurement  of  energy  fluxes;  a  thin  film  heat 
flux  gage,  a  fine  wire  thermocouple  and  a  pyroelectric  sensor.  Each  device  has  its 
own  advantages  and  disadvantages  in  the  measurement  of  laser  energy  flux.  In  this 
appendix,  the  three  energy  flux  measurement  techniques  using  the  three  devices 
will  be  discussed.  The  method  v/hich  employed  a  fine  wire  thermocouple  was  found 
to  be  most  convenient  and  reliable. 

A  thin  film  heat  flux  gage  was  tried  first  among  the  three  devices  to  measure  the 
incident  energy.  Since  the  sensing  area  of  the  thin  film  heat  flux  gage  was  too  large 
to  measure  the  distributions  of  energy  fluxes,  a  mask  with  a  small  hole  (^I.Omm) 
was  used  to  block  the  incident  laser  beam  except  the  portion  which  passes  through 
a  hole  which  has  a  predetermined  diameter  on  the  mask.  An  anodized  aluminum 
mask  block(l"xl"x3/4")  was  mounted  on  a  Medtherm  Rh-Pt  heat  Gux  gage.  Fig. 
Al  shows  the  assembly  of  the  mask  block  and  the  thin  film  heat  flux  gage.  The 
assemly  is  mounted  on  a  positioning  bed  to  scan  the  energy  intensity  distribution. 
Calibration  is  necessary  to  obtain  heat  flux  values  from  the  output  of  the  heat  flux 
gage.  However,  no  absolute  calibration  was  attempted,  since  not  all  the  incident 
laser  is  absorbed  by  the  gage  due  to  nonunity  of  the  surface  absorptivity.  Only 
relative  energy  flux  distribution  was  measured.  One  way  to  estimate  absolute  energy 
flux  distribution  is  through  the  integration  of  the  relative  profile  and  redistribution 
of  the  total  incident  energy.  The  total  energy  of  incident  laser  beam  was  measured 
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by  a  laser  power  probe.  This  value  is  proportionally  distributed  depending  on  the 
local  relative  intensities. 

Another  technique  used  to  measure  the  laser  beam  intensity  profile  was  using 
fine  wire  thermocouple.  A  fine  wire  thermocouple  (4>25fim)  was  installed  as  shown 
in  Fig.  A2.  In  order  to  minimise  the  effects  of  reflected  beam,  a  thin  stainless 
deflector  was  placed  under  the  thermocouple.  A  distance  should  be  kept  between 
the  thermocouple  and  the  deflector  to  prevent  the  effect  of  hot  plume  generated  from 
the  deflector  due  to  laser  heating.  Energy  fluxes  were  deduced  from  the  temperature 
profile  measured  by  the  thermocouple.  As  in  the  thin  film  sensor  technique,  a 
positioning  table  was  used  to  scan  the  beam  profile  with  the  thermocouple.  Figure 
A3  shows  the  intensity  profile  of  the  COj  laser  beam. 

The  last  method  employed  in  this  study  for  measurement  of  incident  laser 
energy  flux  is  using  a  pyroelectric  sensor.  The  schematic  diagram  Fig.  A4  shows 
the  setup  using  a  pyroelectric  sensor  to  measure  beam  intensity  distribution.  A 
PZT(lead  zirconate  titanate)  sensor  wpj  mounted  on  a  support  which  was  attached 
on  a  positioning  table.  The  PZT  sensor  is  directed  to  the  direction  parallel  to  the 
incoming  laser  beam.  The  rise  time  was  10  ns  with  50  ohm  resistor.  The  responsivity 
was  1.5  na/W  and  the  capacitance  of  the  sensor  was  250  pf.  Additional  two  100 
K fl  were  used  to  increase  sensitivity.  The  resistance  was  75KCI  and  the  output 
sensitivity  was  75  mV /W.  A  chopper  is  required  to  measure  intensity  for  a  longer 
duration  than  the  rise  time  which  is  about  100  ns. 
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Fig.  A2  Intensity  Measurement  by  Fine  Wire  Thermocouple 
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X-Y  POSITIONING  BED 


Fig.  A4  Schematic  Diagram  of  Pyroelectric  Sensor  Measurements 
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Uniformity  of  beam  intensity  is  one  of  the  most  important  properties  of  laser 
beam  when  the  laser  beam  is  used  as  a  energy  source  for  ignition  study.  However, 
due  to  the  inherent  structure  of  the  laser,  the  intensity  of  laser  beams  are  generally 
not  uniform.  Unless  the  laser  beam  intensity  is  spatially  uniform,  the  energy  flux 
can  not  be  defined  with  a  single  value. 

Many  different  techniques  were  discussed  to  uniformise  the  intensity  of  laser 
beams  which  were  originally  nonuniform[Al].  Some  of  the  uniformising  techniques 
are  from  solar  energy  research.  A  Kaleidoscope  was  considered  as  the  most  promis¬ 
ing  device. 

The  operation  and  construction  of  kaleidoscopes  were  discussed  extensively 
by  Ohlemiller  and  Summerfield[Al).  In  this  appendix,  only  brief  description  on 
the  principle  of  kaleidoscopes  will  be  presented.  Kaleidoscope  is  a  square  channel 
consisting  of  four  parallel  reflecting  surfaces.  Expanding  laser  beam  which  enters  the 
kaleidoscope  will  be  reflected  off  of  the  inner  walls  at  the  given  angle  of  incidence 
and  continue  through  the  channel.  Since  the  inner  walls  of  the  kaleidoscope  are 
parallel,  the  angle  of  incidence  is  unchanged  as  the  beam  proceeds  through  the 
channel.  Thus,  the  kaleidoscope  breaks  expanding  beam  entering  it  and  recombines 
the  segments  of  the  beam  at  the  exit  of  the  kaleidoscope.  Figure  A5  shows  the 
schematic  diagram  of  the  Kaleidoscope.  The  COj  laser  beam  is  focused  and  diverges 
entering  a  square  mirror  channel.  Some  part  of  the  laser  beam  reflects  from  the 
inner  mirror  surface  two  or  three  times,  while  the  center  portion  of  the  incoming 
beam  does  not  experience  any  reflections  due  to  small  angles  from  the  center  of 
the  beam.  The  beam  emerging  out  of  the  channel,  which  is  composed  of  many 
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different  beam  segments,  is  recombined  by  a  reimaging  lens.  Figure  A5  also  shows 
the  pattern  of  CO2  laser  beam  at  different  locations  after  the  Kaleidoscope.  Due  to 
the  coherence  of  laser  beams,  a  strong  interference  appears  on  the  output  patternc. 
Thus,  it  is  necessary  to  uniformize  the  beams  further  by  using  vibrating  or  rotating 
mirrors  which  average  the  interfered  beam  patterns  timewise. 

Another  device  used  to  achieve  a  uniform  laser  beam  is  a  mask.  As  discussed 
in  Appendix  I,  the  intensity  distribution  of  the  laser  used  in  this  study  is  not  a 
totally  Gaussian.  The  pattern  was  closer  to  TEM10.  The  intensity  at  the  center 
of  the  beam  was  lower  than  those  of  surrounding  locations.  To  invert  this  profile, 
a  stainless  mask(125/im  thick)  with  a  hole(^Smm)  was  placed  on  the  laser  beam 
path.  Due  to  the  converging  effect  of  the  mask,  the  original  intensity  distribution 
became  quite  uniform  in  the  area  of  a  circle  of  radius  4mm  from  the  center.  Figure 
6  shows  the  intensity  distribution  of  the  CO2  laser  beam  uniformized  using  a  mask. 
Figure  A6  shows  the  intensity  distribution  of  the  laser  beam  when  the  beam  mask 
was  used.  Figure  A7  shows  one  cross-sectional  distribution  of  the  laser  intensity. 
Since  the  radius  of  the  sample  was  about  4mm,  the  intensity  of  the  laser  directed 
on  sample  surface  be  ame  uniform  with  a  variation  of  10%. 
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Appendix  C 


ABSORPTION  OF  CO,  LASER  BEAM  ENERGY  BY  CO,  GAS 

The  problems  related  to  the  absorption  of  the  CO2  laser  beam  energy  by  the 
combustion  product  gases  were  addressed  by  Kashiwagi[24]  in  his  ignition  study  of 
PMMA  and  red  oak.  Edwards  et  al.[A2]  studied  the  band  absorption  characteristics 
of  CO2  gas  for  a  range  of  wave  lengths.  The  study  shows  the  existence  of  a  major 
absorption  band  near  10.4  /*m.  As  the  temperature  of  the  gas  increases,  the  width 
of  the  absorption  band  is  widen  so  that  the  absorptance  at  10.6  pm  is  significantly 
increased.  The  absorption  band  of  H2O  near  10.6  /um  shows  a  significant  amount 
of  the  CO2  laser  energy  can  be  absorbed  by  H2O  vapor. 

In  order  to  investigate  the  CO2  laser  absorption  by  CO2  gas,  a  attenua¬ 
tion  test  chamber  was  designed(Fig.  A8).  A  cylindrical  plexiglass  tuking(ID=2", 
OD=2.5",  Invariable)  was  placed  between  two  aluminum  closures.  KC1  windows 
were  mounted  on  the  two  closures,  cne  on  the  top  and  the  other  on  the  bottom  of 
the  chamber.  The  length  of  the  test  chamber  can  be  varied  using  different  plexiglass 
tubing  with  different  length.  Since  the  O-ring  groves  are  on  the  aluminum  closures, 
no  machining  on  the  plexiglass  tubing  is  necessary. 

Figure  A9  shows  the  attenuation  of  the  CO2  laser  beam  intensity  by  C02  gas 
present  in  the  path  of  the  beam.  At  atmospheric  pressure  the  attenuation  increased 
monotonously  between  0%  to  60%  of  CO2  concentration.  As  the  CO2  concentration 
increases  further,  the  slope  of  the  curve  becomes  less  steep.  As  approaching  80%  of 
CO2  gas  concentration,  the  attenuation  reaches  a  limit  value. 

The  attenuation  results  shows  a  significant  amount  of  the  incident  laser  beam 
s  absorbed  by  the  CO2  gas  in  the  beam  path.  One  of  the  main  absorption  bands 
of  CO2  gas  is  at  10.4  nm.  Since  the  wave  length  of  the  CO2  laser  is  at  10.6  /im, 
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a  substantial  portion  of  the  laser  beam  can  be  absorbed  if  the  absorption  band  is 
widened  at  elevated  temperatures.  Water  vapor  is  known  to  be  one  of  the  absorber 
of  the  COa  laser.  However,  absorption  data  for  water  vapor  was  scared.  The 
absorption  test  with  water  vapor  has  not  been  carried  out  in  this  study.  Futher 
studies  are  necessary  for  absorption  of  CO 2  laser  by  other  gases  as  well  as  water 
vapor. 
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Fig.  A8  Beam  Attenuation  Test  Chamber 


APPENDIX  D 


CO,  LASER  OPERATING  PROCEDURE 

1.  Set  the  chiller  temperature  at  15°  C  and  turn  on  the  chiller  power.  Make  sure 
that  the  water  outlet  valve  is  fully  open.  Check  the  water  level  of  the  chiller 
through  the  opening  on  the  top  of  the  chiller.  If  the  water  level  is  lower  than 
1  inch  from  the  top,  refill  it  with  distilled  water.  Do  not  use  the  tap  water. 

2.  Turn  on  the  warning  light  switch  located  on  the  wall  behind  the  CO2  laser. 
Check  if  the  red  warning  light  above  the  entrance  door  is  blinking. 

3.  Unlock  the  breaker  box  (if  locked)  on  the  wall  behind  the  COj  laser  and  turn 
on  the  COa  laser  power  switch. 

4.  Open  the  CO2  laser  mix  bottle,  record  the  pressure  on  the  log  book,  and  set 
the  regulator  to  10  psi. 

5.  Turn  on  the  exhaust  fan  (the  switch  of  the  fan  is  located  on  the  wall  beside  the 
door  of  the  test  cell). 

6.  Turn  on  the  power  of  the  interlocking  system.  If  the  interlocking  system  is  to 
be  defeated,  turn  on  the  interlocking  defeat  switch  on  the  box. 

7.  Make  sure  that  all  the  sliding  glass  windows  on  the  laser  are  properly  closed. 

8.  Check  the  chiller  water  temperature  on  the  dial  gauge  on  the  chiller.  If  the 
temperature  is  higher  than  15°  C,  wait  for  a  few  more  minutes  until  it  reaches 
15°  C. 

9.  Turn  on  the  power  meter  which  is  located  under  the  laser  control  panel. 

10.  Set  the  right  knob  on  the  mode  control  of  the  laser  control  panel  to  CW. 
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11.  Set  the  left  knob  on  the  mode  control  to  OFF. 

12.  Set  the  knob  on  the  current  control  to  halfway  between  0  and  50  mA. 

13.  Set  the  left  keyswitch  to  ON  position. 

14.  Turn  on  the  laser  by  pushing  the  black  pushbutton  switch.  Note  that  the 
orange  laser  lamp  on  the  front  of  the  laser  is  lit. 

15.  Adjust  the  purge  rate  using  the  purge  rate  control  knob  to  maintain  the  tube 
pressure  at  19-20  torr. 

16.  Set  the  left  knob  on  the  mode  control  to  ON. 

17.  Set  the  current  control  knob  to  50  mA  and  check  the  anode  voltage  at  15  KV 
by  depressing  the  tube  drive  switch  to  voltage. 

18.  The  reading  on  the  power  meter  should  be  about  800  Watts.  If  the  power  is 
lower  than  790  Watts,  adjust  the  mirrors  inside  the  laser  following  the  align¬ 
ment  procedure  described  elsewhere  (See  Alignment  Procedure). 

19.  Recycling  mode  can  be  used  to  save  the  laser  mixture  gas.  Follow  the  recycling 
procedure  to  use  the  recycling  mode  (See  Recycling  Procedure). 

20.  Select  the  laser  power  level  by  adjusting  the  current  using  the  current  control 
knob.  Do  not  to  try  to  use  a  power  level  lower  than  100  Watts,  since  at  a  lower 
level  the  laser  becomes  unstable. 

21.  If  a  long  operation  is  required,  use  the  cooling  system  for  the  external  optics 
(See  Water  Cooling  Procedure). 
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Gate  Control  Mode 


22.  The  gate  mode  is  used  to  control  the  laser  generating  time  period.  To  use  the 
gate  mode  control,  set  the  left  knob  on  the  mode  control  to  GATE  position. 

23.  Select  the  desired  laser  generating  time  by  setting  the  gate  control  time  in  the 
combination  of  the  gate  control  knobs.  (For  example,  set  the  left  knob  to  100 
and  the  right  knob  to  4,  then  the  gate  control  time  is  100  x  4  =  400  ms.) 

Continuous  Wave  Operation 

24.  In  order  to  use  the  laser  beam  in  continuous  wave  in  either  the  manual  control 
mode  or  the  gate  control  mode,  set  the  right  knob  on  the  mode  control  to  CW. 
For  the  manual  control,  set  the  left  knob  on  the  mode  control  panel  to  ON 
position.  Now  the  laser  is  on.  To  send  the  laser  beam  out  of  the  laser,  follow 
the  procedure  from  item  26. 

Pulsed  Mode 

25.  To  use  the  pulsed  mode,  set  the  right  knob  on  the  mode  control  to  either  SP 
for  single  pulse  or  REP  for  repeated  pulse. 

26.  Select  the  desired  frequency  and  pulse  length  by  setting  the  frequency  and 
pulse  length  knobs  on  the  pulse  control  panel  in  the  same  way  as  for  the  gate 
control.  For  the  single  pulse,  it  is  not  necessary  to  set  the  knobs. 

Shutter  Control 

27.  To  send  out  the  laser  beam,  set  the  keyswitch  to  the  shutter  power  Opposition. 

28.  Turn  on  the  shutter  power  by  depressing  the  red  button. 
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29.  Set  the  shutter  control  toggle  switch  to  the  OPEN  position.  Now,  the  laser 
beam  is  delivered  out  of  the  laser. 

30.  For  the  gate  control,  turn  on  the  switch  on  the  power  supply  box,  and  turn  it 

off. 

ShutrPpwn  Procedure 

31.  After  the  use  of  the  laser  beam,  set  the  shutter  control  toggle  switch  to  the 
closed  position.  Then,  turn  off  the  keyswitch  for  the  shutter  power. 

32.  Set  the  left  knob  on  the  mode  control  panel  to  the  OFF  position.  Set  the  right 
knob  on  the  mode  control  panel  to  the  CW  position. 

33.  Turn  off  the  laser  by  depressing  the  red  pushbutton  switch  on  the  laser  power 
panel.  Set  the  power  keyswitch  the  OFF  position. 

34.  Turn  off  the  cooling  water  to  the  external  optics  (if  it  was  being  used). 

35.  Turn  off  the  chiller. 

36.  Record  time  and  the  bottle  pressure  in  the  log  book. 

37.  Close  the  main  valve  to  the  gas  bottle  of  the  laser  mix.  Release  the  regulator 
valve. 

38.  Shut  off  the  CO3  laser  power  switch.  Turn  off  the  warning  light  switch. 

39.  Turn  off  the  exhaust  fan  and  the  interlocking  switch. 
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I, _ ,  have  read  and  understood  the  above  operating  procedure 

as  well  as  the  safety  precautions  in  'he  CO2  laser  manual. 


Signature 


Date 


The  above  person  is  authorised  to  operate  the  CO2  laser. 


Prof.  K.  K.  Kuo 


Date 


cc:  Mr.  John  Raiser,  Mech.  Engr. 

Mr.  Robert  Houtz,  College  of  Engr. 

Mr.  H.  O.  Triebold,  Jr.,  Mgr.,  Safety  Division. 
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APPENDIX  E 


SPECIFICATIONS  OF  EQUIPMENT  USED  IN  LASER  IGNITION  -STUDY 

•  Laser 

•  Model  Super48  COj  Laser,  800  Watt  in  continuous  wave,  3500  Watt  in  pulsed 
mode. 

.  Coherent,  3210  Porter  Drive,  Palo  Alto,  CA  04304,  (415)  493-2111. 

•  Thermal  Image  Plate 

-  Thermal  Image  Plate  and  Ultraviolet  Lamp. 

.  Optical  Engineering,  Inc.,  P.O.Box  696,  Coffey  Lane,  Santa  Rosa,  CA 
95402,  (707)  528-1060. 

•  Laser  Power  Probe 

-  Model  P-20Y  and  P-1000Y  Laser  Power  Probe. 

.  Optical  Engineering,  Inc.,  P.O.Box  696,  Coffey  Lane,  Santa  Rosa,  CA 
95402,  (707)  528-1060. 

•  Mirrors  and  Lenses 

-  3"  Silicon  Mirrors  with  Enhanced  Silver  Coating. 

-  1.1"  dia.  Plano-Convex  Lenses  with  AR  Coating. 

.  Two-Six  Incorporated,  Saxonburg  Boulevard,  Saxonburg,  PA  16056,  (412) 
352-1504. 

.  Laser  Optics,  Inc.,  P.O.Box  127,  Danbury,  Conneticut  06810,  (203)  744- 
4160. 

•  Windows 

-  3"  Zinc  Selenide(ZnSe)  Window  with  AR  Coating. 

.  Two-Six  Incorporated,  Saxonburg  Boulevard,  Saxonburg,  PA  16056,  (412) 
352-1504. 

-  1/2"  KC1  Window  for  Low  Pressure  Chamber. 

.  Harshaw  Chemical  Co.,  6801  Cochran  Road,  Solon,  Ohio,  44139,  (216) 
248-7400. 

•  Optic  Mounts 

-  Model  LP-1  Five  Axis  Gimbal  Optic  Mounts 

.  Newport  Corporation,  18235  Mt.  Baldy  Circle,  Fountain  Valley,  CA 
92708,  (714)  963-9811. 

-  Dovetail  Optical  Track  and  Slides. 

.  Thomas  Tool  and  Die,  Inc.,  8805  Bradley  Avenue,  Sun  Valley,  CA  91352, 
(818)  767-6511. 
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U nitron  ZST  with  Fiber  Optic  Ring  Illuminator. 

.  Unitron  Inc.,  175  Express  Street,  Plainview,  N.Y.,  11803,  (516)  822-4601. 

Ehatadtoato 

Model  SGD-040L  Fast  Response  Silicon  Photodiode. 

.  EG&G  Photon  Devices,  35  Congress  Street,  P.O.Box  5006,  Salem,  MA 
01970-6526,  (617)  745-3200. 

Pyroelectric  Detector 

Lead  Zirconate  Titanate  (PZT)  Laser  Pvlse  Detector. 

.  Barnes  Engineering  Co.,  30  Commerce  Road,  Stamford,  Conneticut  06904, 
(203)  348-5381. 


25 fxm  Pt-Ptl3%Rh  Fine  Wire  Thermocouple. 

.  Omega  Engineering,  Inc.,  P.O.Box  2669,  Stamford,  CT  06906,  (203)  359- 
1660. 

5 fim  Pt-Ftl0%Rh  Wollaston  Process  Wires. 

.  Sigmund  Cchn  Corp.,  121  S.  Columbus  Avenue,  Mount  Vernon,  NY  10553, 
(214)  664-5300. 

Model  TCF30?-tth-Pt-M  Thin  Film  Thermocouple. 

.  Medtherm  Corporation,  P.O.Box  412,  Huntsvill  Alabamb,  35804,  (205) 
837-2000. 


Digital  Oscilloscope 

Nicolet  Series  2090  Oscilloscope. 

.  Nicolet  Instrument  Corporation,  5225-2  Verona  Road,  Madison,  Wisconsin 
53711.  (608)  271-3883. 

High  Speed  Movie  Camera 

HYCAM  16mm  High  Speed  Movie  Camera,  Model  41-0006. 

.  Redlake  Corr.,  15005  Concord  Circle,  Morgan  Hill,  CA  95037,  (408)  779- 
6464. 


High  Speed  Video  Camera 

Spin  Physics  High-Speed  Motion  Analyzer. 

.  Spin  Physics  Division,  Eastman  Kodak  Co.,  3099  Science  Park  Rd.,  San 
Diego,  CA  92121,  (619)  481-8182. 

Videoeranh  Printer 

AIDI  Videograph  Printer  CT1500. 

.  Advanced  Imaging  Devices,  Inc.,  875  Mandc  Avenue,  Mountain  View,  CA 
94943. 


Video  System 

Panasonic  PK-956  Video  Camera 
Panasonic  PV-5500  Video  Cassette  Recorder 

.  Panasonic  Co.,  One  Panasonic  Way,  Secaucus,  New  Jersey  07094,  (201) 
348-7000. 


IBM-PC. 

.  International  Business  Machines  Corporation,  P.O.Box  1328-S,  Boca  Ra¬ 
ton,  FL  33432. 

Schlieren  System 

Light  Source(100W  Tungsten-Halogen  Lamp  28-8472) 

Condensor  Lens(D=18mm,  f=40rcm  Achiomat) 

Source  Spatial  Filter(Rectangular  Aperture  23-1928} 

Focussing  Lens(D=40mm,  f=90mm  Achromat) 

.  Ealing  Electro-Optics  Inc.,  22  Pleasant  St.,  South  Natict,  MA  01760, 
(617)651-8100. 

Zeiss-Type  Optica!  V-Bench(F83-008) 

Transverse  and  Vertical  Slides(22-4899)  and  Base(22-4170) 

60mm  Pin  Carriers  (F60-796) 

Mounting  Pins(F41-229) 

Screen  Holder(22-8700) 

Rectagular  Ground  Glass(F42-785) 

.  Edmund  Scientific  Co.,  101  E.  Gloucester  Pike,  Barrington,  N.J.  08007, 
(609)  547-3488. 

Gas  Control  System 

High  Pressure  Tubings,  Valves,  and  Connectors. 

.  High  Pressure  Equipment  Company,  Inc.,  1222  Linden  Avenue,  Erie,  PA 
16505,  (814)  838-2028. 

Coroprewor 

Single-End  Compressor. 

.  Newport  Scientific,  Inc.,  8246-E  Sandy  Court,  Jessup,  MD  20794-0189, 
(301)  498-6700. 

Vacuum  Pump 

Model  KC-3  Kinney  Vacuum  Pump. 

.  Kinney  Vacuum  Company,  495  Turnpike  Street,  Canton,  MA  02021,  (617) 
828-9500. 


•  Pressure  Regulator 

•  Matheson  Model  3060  Series  High  Pressure  Regulators. 

.  Matheson  Gas  Products,  P.O.Box  1S87,  Secaucus,  N.J.  07094. 

•  Gas  Chromatograph 

-  Varian  3700  Gas  Chromatograph. 

.  Varian  Instrument  Group,  25  Hanover  Road,  Florham  Park,  N.J.  07932, 
(201)822-3701. 

-  Columns  and  Syringes. 

.  Supelco  Park,  PA  16823-0048,  (814)  359-3441. 

-  Standard  Gases. 

.  Scott  Specialty  Gases,  Route  611,  Plumsteadville,  PA  18949,  (215)  766- 
8861. 

.  Union  Carbide  Corporation,  Linde  Division,  National  Specialty  Gases  Of¬ 
fice,  100  Davidson  Avenue,  Somerset,  NJ  08873,  (201)  271-2600. 

•  Strobe  Light 

-  Strobotac  Model  1531. 

.  General  Radio  Co.,  300  Baker  Avenue,  Concord,  MA  01742,  (617)  646- 
0550. 
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APPENDIX  F 


TEST  PROCEDURES  AND  CHECKLIST  FOR  CRACK  PROPAGATION  TESTS 

A.  Igniter  Assembly: 

_  1.  Cut  igniter  charge  and  record  mass  and  propellant  type. 

_  2.  Select  a  primer  with  a  resistance  20  -  SO  ohms  and  record. 

_  3.  Assemble  primer,  primer  holder  and  lexane  insulator. 

_  4.  Check  primer  resistance.  If  it  has  changed  significantly, 

check  leads  and  reassemble. 

_ 5.  Place  desired  amount,  of  propellant  shavings  in  and  around 

the  primer. 

_  6.  Assemble  remaining  igniter  pieces  with  the  yellow  marks 

alligned. 

_  7.  Check  igniter  circuit  continuity. 

_  8.  Short  igniter  leads. 

B.  Sample  preparation: 

_  1.  Measure  and  record  the  mass  and  dimensions  of  the  sample. 

_  2.  Paint  a  dot  on  the  front  surface  of  the  propellant  for 

focusing  purposes. 

_  3.  If  Spin  Physics  camera  will  be  used,  mark  the  crack  tip 

location  with  paint  and  measure  distance  between  crack  tip  and 
focusing  dot. 

C.  Chamber  Asseraly: 

_  1.  Install  the  desired  Ng  inlet  lines  (except  on  igniter)  and 

check  for  unobstructed  flow. 

_  2.  Assure  sacrificial  and  main  windows  properly  fit  into  test 

chamber  and  clean  window  surfaces. 


_  3.  Spread  flame  retardant  grease  over  all  sample  surfaces 

(except  or ack  surface)  and  Install  in  test  chamber. 

_  4.  Insure  propellant  sample  will  not  extend  beyond  sample 

holders  when  compressed. 

_ _  5.  Install  main  O-ring  with  joined  ends  at  the  rear  of  the 

chamber.  Apply  permatex  to  O-ring  and  groove  * 

_  6.  Install  sacrificial  window  against  the  propellant  sample. 

_  7.  Install  the  main  window  into  the  chamber  window  cover. 

_  8.  Bolt  the  chamber  window  assembly  to  the  main  chamber, 

being  careful  to  evenly  compress  the  sample. 

_  9.  Measure  and  record  top  bursting  diaphragm  thickness  and 

install  diaphragm  and  hood  assembly  to  top  of  chamber.  Check  for 
proper  allignment  of  diaphragm  holder  to  allow  maximum  flow  area. 

- -  10.  Check  side  rupture  disk  and  replace  if  necessary. 

_  11.  Install  igniter  assembly  to  bottom  of  chamber  being  sure 

to  alllgn  the  yellow  marks.  Check  igniter  resistance  and  short  the 
igniter  circuit. 

_  12.  Connect  Ng  line  to  igniter.  If  not  used,  ensure  that 

igniter  ports  are  plugged. 

D.  Instrumentation: 

_  1.  Verify  that  all  igintion  circuits  are  in  the  off  position. 

_  2.  Install  pressure  transducers  and  record  their  respective 

locations,  sensitivities,  and  data  lines. 

_  3.  Connect  water  lines  and  check  for  leaks.  Turn  off  water 

supply. 

_  4.  Set  up  Hycam  and  install  dummy  film.  Insure  the  film 

speed  is  reduced  for  ignition  circuit  test. 


_  5.  Connect  volt-meter  to  ignition  box  and  run  ignition  test. 

_  6.  Check  focus  of  the  Hycam. 

_ _  7.  Load  new  film,  depress  microswiteh,  and  insure  triggering 

arm  operates  properly  as  the  camera  door  is  secured.  Reset  PPS 
dial. 

_  8.  Check  and  record  all  camera  settings. 

_  9.  Record  test  cell  temperature. 

_  10.  Setup  and  focus  Spin  Physics  camera  (if  used). 

_  11.  Remove  lights  used  for  focusing. 

_  12.  Install  exaust  fan  and  connect  to  Lights  No.  1  receptacle. 

_  13.  Check  data  line  connections  to  charge  amplifiers. 

_  14.  Set  input  voltages  and  calibrate  charge  amplifiers. 

_  15.  Connect  charge  amplifier  output  to  Physical  Data  recorder 

and  record  settings.  Turn  on  Physical  Data  unit  and  hit  reset 
button. 

E.  Test  firing  procedure: 

1.  Perform  final  check  of  all  connections. 

_ Hycam  _ Spin  Physics 

_  Transducers  _  Water  connections 

_ Exhaust  fan  _ Ng  lines 

_ Physical  Data  _ Charge  Amps 

_  2.  Open  outside  door. 

_  3.  Secure  safety  key. 

_  4.  Connect  igniter  leads  to  igniter  box. 

_  5.  Check  igniter  resistance. 

_  6.  Close  and  secure  test  cell  door. 

_  7.  Turn  on  warning  light  and  sound  test  alarm. 
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.  8.  Verify  eharge  amps  are  in  CHARGE  position. 

.  9.  Turn  on  igniter  power  and  charge  capacitor. 

10.  Turn  on  ignition  panel  switches  3  thru  6  and  event  switch. 

11.  Turn  off  safety  switch. 

12.  Unground  charge  amplifiers,  reset  Physical  Data  recorder. 

13.  Turn  on  Spin  Physics  Camera  and  wait  for  ready  light. 

14.  Turn  on  ignition  switch. 

15.  After  rupture  disk  bursts,  inject  Ng. 

16.  When  propellant  has  stopped  burning,  stop  Spin  Physics 


camera. 

F.  Post  Firing  Procedures: 

_  2.  Ground  charge  amplifiers. 

_  3.  Turn  off  all  ignition  switches  and  turn  on  exhaust  fan 

switch. 

_ 4.  Make  a  hard  copy  of  pressure  traces  on  X-v  il,-'-.r,er. 

_  5.  Transfer  pressure  data  to  Nicolet  end  store  on  floppy 

disk. 

_  6.  Generate  calibration  signals  for  transducers  and  record  on 


X-Y  plotter 


DATA  SHEET  FOR  CRACK  PROPAGATION  EXPERIMENTS 


CRACK  PROPAGATION  TEST  NO. 

DATE 

participants: 

9  i 

i  > 

TEMPERATURE:  .  T(room) 

C 

T(cell)  C 

CRACK  GEOMETRY  AND  SAMPLE  SPEC.: 

IGNITER: 

Propellant  type 

Primer 

Configuration  _ 

Propellant  type 

Length 

• 

Dimensions 

Kax  gap  width 

Mass 

Inner  taper  angle 

Resistance 

Taper  heigth 

Thickness 

FLYING  PIN  IGNITER: 

Height 

Primer  resistance 

Fuse  Hire  Location 

, 

Propellant  weight 

Propellant  type 

CAMERA: 

Diaphragn  thickness 

Type  of  Lens 

Diaphragm  diameter 

Framing  Rate 

Frames/s 

Diaphragm  material 

Optical  Head:  1/A 

,1/2  ,  frame 

PPS  Dial 

RUPTURE  DISK  ACTUATOR: 

Multiplier 

Primer  resistance 

Time  Marker 

Pulses/s 

Propellant  weight 

Remote  firing  at 

ft 

Propellant  type 

Distance  of  focu.> 

ft 

Brake  Dial 

FILM: 

Servo 

ASA  NO. 

H/L  speed  reducer 

Length 

F.  Stop: 

Color  ,  B  &  W 

APPENDIX  G 


ERROR  ANALYSIS 

In  the  multiple  regression  of  the  data  in  Table  VI,  the 
exponential  relationships  assumed  for  the  dimensionless  parameters 
[equations  (3-36)  through  (3-39)]  were  rearranged  into  a  linear 
relationship  based  on  the  natural  logarithm  of  each  value.  Using  the 
basic  properties  of  logrithmic  representations,  the  uncertainty  in  each 
dimensionless  quantity  can  be  expressed  as 

ln(xi)  a  ln(x)  +  ex/x  (G-41 ) 

where  Xj  is  the  actual  value,  x  is  the  measured  value  and  ex  is  the 

81 

error  interval  associated  with  that  measurement.  The  method  of  least 
square  analysis  used  on  the  measured  data  will  give  more  accurate 
results  when 

ln(x)  »  ex/x  (G-42) 

An  expression  for  the  relative  error  in  each  dimensionless  parameter  is 

i*x  =  ex/(xln(x))  (G-43) 

Table  VII  contains  estimations  of  the  various  errors  involved  in 
measuring  each  of  the  variables  in  the  experimental  setup.  Equation 
(G-43)  can  then  be  evaluated  for  each  of  the  dimensionless  parameters. 
This  equation  indicates  that  for  values  of  x  which  approach  unity,  the 
relative  error  could  be  significant.  Since  the  error  of  c  (n)  is 
assumed  to  be  zero,  its  correlation  is  not  affected.  However,  a^  does 
have  a  finite  error  associated  with  it.  For  the  worst  case  of 
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■T  «  1.1,  rx  *  0.095  or  9.5%  relative  error.  This  could  be  Improved  by 
referencing  a^  to  •  temperature  outside  the  normal  operating  range. 

Using  equation  (G-43)  and  the  values  listed  in  Table  4,  the 
relative  error  was  calculated  for  eaoh  dimensionless  parameter  based  on 
its  average  value. 


A 


Vo 


(G-44) 


Assume  e.  «  e.  therefore  e  /  .  *  e(.  )/Ln.  Based  on  the  average 
LD  LD  A  A  LD  0 

value,  .  «  0.326,  r  *  9.91. 

*  A 

Using  this  method,  the  following  values  of  rx  are  obtained: 


9.9% 

3.9% 

0.0% 

0.3% 

1.7% 

3.4% 

=  5.9% 
=  8.2% 
2.1% 
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Table  VII 


Estination  of  Measurement  Errors 


Physical  Error 

Variable  Souroe  and  Type  of  Error _ Interval  (♦/-) 


1 

0 

Random  error  In  the  measurement  and 
calculation  of  the  Inner  taper  angle. 

o 

rw 

• 

o 

The  systematic  error  of  symmetry 

0.5° 

•t 

Random  error  In  the  measurement  of  the 
room  temperature 

0.5  K 

Systematic  error  in  the  difference  between 
romm  temp,  and  initial  propellant  temp. 

3.0  K 

Random  and  systematic  error  in  converting 
temperature  to  a.  (note:  in  the  temperature 
range  of  interest  da^/dT  is  near  zero) 

0.01 

Accuracy  of  the  a_  data  (note:  a.  is 
defined  near  operating  temperatures) 

NA 

n 

Since  n  is  a  discrete  quantity*  it  is  assumed 
to  have  no  random  error 

0.0 

•o 

Random  error  of  the  measuring  instrument 

0.25  nan 

h> 

Random  error  in  determining  the  bright  flame 
front  tip  location 

2.0  mm 

Systematic  error  of  correlation  between  the 
flame  front  and  actual  tip  location 
(note:  assumed  to  be  small) 

unknown 

V 

P 

Random  error  of  measuring  the  change  of 
flame  front  location  with  respect  to  time 
(due  to  the  averaging  effect) 

0.0 

P 

Systematic  error  of  the  installation,  wear, 
calibration  and  recording  of  pressure  data 

5S 

dp/dt 

Since  dp/dt  is  calculated  using  a  very  large 
dp  and  t  is  very  accurate  the  error  should  not 
be  greater  than  in  P 

5* 
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